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ABSTRACT 


Multispectral  imaging  (MSI)  data  collected  at  multiple  angles  over  shallow  water  provide 
analysts  with  a  unique  perspective  of  bathymetry  in  coastal  areas.  Observations  taken  by 
DigitalGlobe’s  WorldView-2  (WV-2)  sensor  acquired  at  39  different  view  angles  on  30 
July  2011  were  used  to  determine  the  effect  of  acquisition  angle  on  derived  depth.  The 
site  used  for  this  study  was  on  the  island  of  Oahu,  focused  on  Kailua  Bay  (on  the 
windward  side  of  the  island).  Satellite  azimuth  and  elevation  for  these  data  ranged  from 
18.8  to  185.8  degrees  and  24.9  (forward-looking)  to  24.5  (backward-looking)  degrees 
(respectively)  with  90  degrees  representing  a  nadir  view.  Bathymetry  were  derived 
directly  from  the  WV-2  radiance  data  using  a  band  ratio  approach.  Comparison  of  results 
to  LiDAR-derived  bathymetry  showed  that  varying  view  angle  impact  the  quality  of  the 
inferred  bathymetry.  Derived  and  reference  bathymetry  have  a  higher  correlation  as 
images  are  acquired  closer  to  nadir.  The  band  combination  utilized  for  depth  derivation 
also  has  an  effect  on  derived  bathymetry.  Four  band  combinations  were  compared,  and 
the  Blue  &  Green  combination  provided  the  best  results. 
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I.  INTRODUCTION 


A.  PURPOSE  OF  RESEARCH 

Imagery  data  acquired  from  satellites  are  widely  useful  in  the  field  of 
oceanography.  These  data  contribute  to  a  better  understanding  of  the  dynamics  of  ocean 
circulation,  are  useful  for  monitoring  climate  change,  can  be  used  for  navigation  and 
fisheries  management,  and  are  also  helpful  for  improving  models  of  ocean  circulation, 
air-sea  interaction,  weather  forecasting,  and  climate  [ Sanford  et  a/.,  2011]. 

Information  about  shallow  water  bathymetry  is  beneficial  to  scientists  or  groups 
that  require  knowledge  of  ocean  depths  in  a  particular  coastal  location.  The  use  of 
multispectral  imagery  (MSI)  data  has  been  shown  to  adequately  determine  depths  of 
remote  coastal  areas,  as  SOund  Navigation  And  Ranging  (SONAR)  or  Light  Detection 
And  Ranging  (LiDAR)  soundings  or  other  bathymetric  data  may  not  be  available  prior  to 
arrival. 

There  are  many  benefits  to  using  remote  sensing  data  acquired  from  satellites. 
Spaceborne  sensors  can  collect  spectral  data  over  extremely  large  areas  which  can  be 
advantageous  to  regions  not  accessible  on  foot,  or  for  areas  denied  by  hostile  forces. 
These  data  can  also  be  acquired  much  more  rapidly  due  to  the  number  of  space  platforms 
collecting  on  a  regular  basis.  Finally,  improvements  in  collection  capabilities  steadily 
grow  as  space  platform  production  increases  and  sensor  technologies  evolve.  This 
capability  is  extremely  valuable  to  the  Navy  and  Marine  Corps,  who  routinely  conduct 
and  plan  worldwide  amphibious  operations.  Coastlines  are  highly  dynamic  environments 
and  mission  success  is  dependent  upon  knowing  the  underwater  terrain.  Being  able  to 
acquire  bathymetry  in  a  contested  or  denied  region  greatly  improves  the  probability  of 
success  for  these  amphibious  operations. 

The  purpose  of  this  research  was  to  integrate  the  use  of  multiple  satellite  image 
acquisition  angles  over  one  location,  and  detennine  what  role  these  varying  angles  play  in 
bathymetric  depth  determination.  Analyses  of  39  WorldView-2  (WV-2)  images  acquired 
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over  Kailua  Bay  on  the  windward  side  of  Oahu,  Hawaii  were  used  to  reach  a  conclusion 
about  the  accuracy  of  bathymetric  derivation  from  MSI. 


B.  SPECIFIC  OBJECTIVES 

The  objective  of  this  study  was  to  test  the  potential  of  bathymetric  derivation 
using  WV-2  imagery  acquired  at  multiple  angles,  and  then  report  upon  the  role  that 
image  acquisition  angle  plays  in  depth  detennination. 

The  motivation  for  this  work  originated  from  the  need  to  determine  bathymetry 
from  only  a  single  spectral  image  of  the  coastal  region  in  question.  There  is  no  guarantee 
that  this  image  will  have  been  acquired  at  optimal  viewing  geometry,  e.g.,  nadir.  Without 
sufficient  time  to  task  a  satellite  and  acquire  data  at  a  particular  viewing  angle  before  a 
site  visit,  results  from  this  research  strive  to  provide  a  better  understanding  of  how  to 
manipulate  image  data  in  order  to  obtain  a  better  understanding  of  bathymetry.  Analysis 
of  multi-angle  MSI  data  was  used  to  quantify  the  effects  of  varying  satellite  acquisition 
angle  on  accurate  determination  of  bathymetry. 
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II.  BACKGROUND 


This  research  focused  on  the  use  of  passive,  optical  remote  sensing  systems  for 
image  acquisition.  Passive  sensors  are  those  that  rely  on  incoming  solar  radiation  to 
illuminate  the  targets  on  land  and  in  the  water  [ Camacho ,  2006].  Optical  sensors  are 
those  that  focus  on  the  visible  (VIS),  near  infrared  (NIR),  and  shortwave  infrared  (SWIR) 
portions  of  the  electromagnetic  (EM)  spectrum  to  observe  radiation  from  targets 
[Camacho,  2006].  The  following  sections  will  explain  these  principles  in  more  detail,  as 
well  as  further  discuss  some  issues  that  affect  remote  sensing  of  ocean  environments,  as 
well  as  imagery  acquired  at  multiple  angles. 

A.  PRINCIPLES  OF  RADIATIVE  TRANSFER 

A  remote  sensing  instrument  or  sensor  receives  energy  that  is  reflected  from  the 
surface  of  the  Earth.  This  energy  is  affected  by  interactions  of  light  with  the  atmosphere 
and  water,  as  well  as  any  particulate  matter  in  the  water  column  [ Camacho ,  2006].  It  is 
imperative  to  understand  the  basics  about  the  EM  spectrum,  radiation,  the  variety  of 
possible  interactions,  and  line  spectra. 

1.  Electromagnetic  Spectrum 

The  EM  spectrum  is  a  conglomeration  of  a  number  of  classifiable  spectral 
regions.  These  include:  gamma  rays,  X-rays,  ultraviolet  (UV)  light,  VIS  light,  IR  light, 
microwaves,  and  radio  waves  (Figure  1)  [Olsen,  2007]. 
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Figure  1.  The  Electromagnetic  Spectrum  (From  Jensen  [2007]) 

This  research  focuses  on  the  VIS  and  NIR  portions  of  the  spectrum,  extending 
from  approximately  400  nm  to  750  nm,  and  750  nm  to  1000  nm  (respectively). 

2.  Spectral  Signatures 

Single  atoms  or  molecules  emit  light  in  the  form  of  line  spectra,  and  an  atom  that 
is  well  isolated  will  radiate  a  discrete  set  of  frequencies  called  a  line  spectrum.  The 
wavelengths  radiated  and  absorbed  are  specific  to  that  atom  or  molecule,  and  are 
presented  as  spectral  signatures.  These  signatures  are  used  to  determine  the  composition 
of  radiating  or  absorbing  gases  and  other  materials  [Olsen,  2007]. 
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Spectral  features,  such  as  those  shown  in  Figures  2  and  3,  enable  analysts  to  better 
differentiate  between  materials  on  the  ground.  Based  on  specific  values,  peaks,  and 
troughs  within  the  graph,  one  can  determine  the  basic  makeup  of  a  material  within  the 
image. 


A  (%) 


Figure  2.  Hyperspectral  Imagery  (HSI)  provides  more  complete  spectral  signatures  than 
MSI  because  there  are  a  much  larger  number  of  bands.  Albedo  values  for 
different  forms  of  algae  and  coral  are  shown  (From  Maritorena  et  al.  [1994]; 

Camacho  [2006]) 
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Figure  3.  MSI  has  spectral  signatures  that  look  less  complete  compared  to  HSI  because 
there  are  fewer  bands  (spectra  are  taken  from  a  scene  used  in  this  research  and  are 

in  nanometers) 


3.  The  Four  Fundamental  Energy  Interactions  with  Matter 

Electromagnetic  radiation  may  be  transmitted,  reflected,  scattered,  or  absorbed. 
The  proportions  to  which  these  interactions  occur  depend  on  the  compositional  and 
physical  properties  of  the  medium,  the  wavelength  or  frequency  of  the  incident  radiation, 
and  the  angle  at  which  the  incident  radiation  strikes  a  surface  [Avery  and  Berlin,  1992; 
Olsen,  2007]. 


a.  Transmission 

With  transmission,  shown  in  Figure  4,  incident  radiation  passes  through 
matter  without  measurable  attenuation.  Different  densities  of  the  material,  however,  can 
cause  radiation  to  be  refracted  or  deflected  from  a  straight-line  path,  and  will  also  alter 
the  velocity  and  wavelength.  The  change  in  EM  radiation  velocity  can  be  further 


6 


explained  by  the  index  of  refraction,  which  is  the  ratio  between  the  velocity  of  the  EM 
radiation  in  a  vacuum  (a  perfectly  transparent  medium)  and  its  velocity  in  a  material 
medium. 


The  index  of  refraction  for  a  vacuum  is  equal  to  1 ,  and  can  never  be  less 
than  1  for  any  substance  [ Avery  and  Berlin,  1992]. 


Figure  4.  Transmission  (After  Avery  and  Berlin  [  1 992]) 

b.  Reflection 

Reflection,  or  specular  reflection,  occurs  when  surfaces  are  smooth 
relative  to  the  wavelengths  of  incident  radiation  (Figure  5).  Specular  reflection,  the 
process  by  which  incident  radiation  bounces  off  the  surface  of  a  substance  in  a  single  and 
predictable  direction,  causes  no  change  to  either  the  EM  radiation  velocity  or  wavelength 
[Avery  and  Berlin,  1992]. 
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Figure  5 .  Reflection  (After  Aveiy  and  Berlin  [  1 992]) 

c.  Scattering 

Scattering,  or  diffuse  reflection,  takes  place  when  incident  radiation  is 
dispersed  in  unpredictable  directions.  Scattering  occurs  when  surfaces  are  rougher 
relative  to  the  wavelengths  of  incident  radiation,  as  shown  in  Figure  6.  The  velocity  and 
wavelength  of  EM  waves,  however,  are  not  affected  by  scattering  [Avery  and  Berlin, 
1992]. 


Figure  6.  Scattering  (After  Avery  and  Berlin  [  1 992]) 
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(1)  Bidirectional  Reflectance  Distribution  Function  (BRDF).  The 
Bidirectional  Reflectance  Distribution  Function  (BRDF)  describes  the  scattering 
characteristics  of  a  material  by  describing  how  light  from  a  source,  incident  on  a  target,  is 
reflected  in  a  given  direction.  In  order  to  calculate  BRDF,  one  looks  at  the  ratio  of  the 
energy  scattered  by  the  target  in  a  particular  direction,  dependent  on  wavelength,  over  the 
energy  that  was  incident  on  the  target  from  a  particular  direction.  All  combinations  of 
inbound  and  outbound  energy  directions  are  then  integrated.  BRDF  is  dependent  on  both 
the  angle  of  incidence  and  on  the  angle  of  reflectance.  Angles  are  defined  in  a  Cartesian 
coordinate  system  by  a  polar  angle,  9,  measured  from  the  surface  normal,  as  well  as  an 
azimuthal  angle,  (p,  measured  from  the  x-axis  (Figure  7)  [ McConnon ,  2010]. 


y 


Figure  7.  Reflection  geometry  used  in  the  definition  of  BRDF 

(From  McConnon  [2010]) 


d.  Absorption 

Absorption  is  the  process  by  which  incident  radiation  is  taken  in  by  a 
medium.  This  occurs  when  a  substance  is  opaque  to  the  incident  radiation,  a  portion  of 
which  is  converted  to  internal  heat  energy,  then  emitted  or  reradiated  at  longer  thermal 
infrared  wavelengths  (Figure  8)  [Avery  and  Berlin,  1992]. 


9 


Emission 


Figure  8.  Absorption  (After  Avery  and  Berlin  [1992]) 

B.  INTERACTIONS  OF  LIGHT  WITH  THE  ATMOSPHERE 

There  are  three  limiting  factors  that  the  atmosphere  introduces  into  the  field  of 
remote  sensing.  These  include  atmospheric  absorption,  scattering,  and  turbulence  [Olsen, 
2007].  Atmospheric  turbulence,  however,  will  not  be  discussed  here  because  its  impact  is 
greater  for  telescopes  looking  up  from  the  earth  through  the  atmosphere  than  for  sensors 
looking  down.  Also,  it  is  important  to  keep  in  mind  that,  overall,  the  atmosphere  is  more 
transparent  in  the  long-wave  IR  (1100  to  1200  nm)  than  in  the  VIS  spectrum  (400  to  750 
nm)  [Olsen,  2007]. 

1.  Atmospheric  Absorption 

Atmospheric  absorption  is  mainly  dependent  on  wavelength,  and  is  most  affected 
by  water  (between  500  and  700  nm),  carbon  dioxide,  and  ozone  (near  1000  nm)  [Olsen, 
2007].  Even  after  sensor  calibration  has  occurred,  the  shape  of  the  measured  spectra 
relative  to  the  underlying  reflectance  spectra  are  altered  because  of  atmospheric 
absorption  lines  [Eismann,  2012].  Common  absorption  spectra  can  be  seen  in  Figure  9. 
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Figure  9.  The  first  four  graphs  show  the  absorption  characteristics  of  N2O, 

02  and  03,  C02,  and  H20,  while  the  bottom  graphic  depicts  the  cumulative  result 
of  all  these  constituents  being  in  the  atmosphere  at  one  time  (From  Jensen  [2007]) 

2.  Atmospheric  Scattering 

Atmospheric  scattering  is  primarily  caused  by  collisions  between  aerosols  and 
particulates  like  dust,  fog,  and  smoke  in  the  atmosphere  [Olsen,  2007].  There  are  two 
main  consequences  of  atmospheric  scattering:  reduction  of  radiant  energy  and  unwanted 
gain  in  the  sensor  [Martin,  2004;  Camacho,  2006].  Scattering  can  also  be  divided  into 
two  types:  Rayleigh  and  Mie  scattering.  Rayleigh,  or  molecular,  scattering  is  primarily 
caused  by  oxygen  and  nitrogen  molecules,  i.e.,  particles  whose  effective  diameters  are 
less  than  the  wavelengths  of  interest.  As  the  size  of  the  particle  increases,  the  scattering 
processes  move  towards  Mie  scattering.  This  type  of  scattering  occurs  when  the  diameter 
of  particles  is  similar  to  the  wavelengths  of  the  energy  being  radiated  [Camacho,  2006; 
Olsen,  2007]. 
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3.  Atmospheric  Compensation  for  Spectral  Imagery 

The  effects  the  atmosphere  has  on  incoming  radiation  reaching  a  remote  sensing 
sensor  results  from  a  variety  of  processes.  Five  possibilities  are  shown  in  Figure  10  [Kay 
et  al.,  2009]. 


— ►  A  Single  or  multiple  scattering  in  the  atmosphere,  by  molecules  or  aerosols 

— ►  B  Scattering  from  the  atmosphere  to  the  water  surface  and  then  reflection  to  the  detector  often 
termed  “sky  glint” 

— ►  C  Reflection  from  whitecaps  on  the  sea  surface 

— ►  T)  Specular  reflection  from  water  surface,  with  direct  transmission  through  the  atmosphere  from 
the  sun  to  the  surface  and  from  the  surface  to  the  detector -  this  is  termed  sun  glint'  in  the 
context  of  this  review. 

— ►  F.  Transmission  through  the  atmosphere  and  air-water  interface  followed  by  scattering  or 
reflection  below  the  water  surface  and  transmission  back  through  the  atmosphere  to  the 
detector. 


Figure  10.  Diagram  showing  routes  by  which  light  can  reach  a  remote 

sensing  detector  (From  Kay  et  al.  [2009]) 


Kay  et  al.  [2009]  assumes  that  if  these  five  processes  are  predominantly 
responsible  for  the  sensor-received  signal,  then: 


I  atm  "t"  nsky  +  TL 


^whitecap 


+  II  glint  +  TL 


water 


(i) 


where  T  is  the  transmittance  of  the  atmosphere  along  the  sensor  view  direction.  The  path 
I  atm  is  the  radiance  that  arrives  at  the  sensor  via  atmospheric  paths  (single  or  multiple 
scattering  in  the  atmosphere,  by  molecules  or  aerosols).  Other  terms  Is ky,  Iwhitecap,  Igiint, 
and  Iwater  are  the  radiances  just  above  the  water  surface  for  light  traveling  by  sky  glint, 
whitecap,  sun  glint,  and  water-leaving  routes.  These  terms  depend  on  the  wavelength, 
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along  with  other  factors.  The  term  Lwa ter  contains  the  infonnation  about  water  column 
and  benthic  features,  and  needs  to  be  separated  from  the  other  terms  if  this  infonnation  is 
going  to  be  retrieved  [Kay  et  al.,  2009]. 

C.  INTERACTIONS  OF  LIGHT  AND  WATER 

Optical  oceanography  is  vital  for  addressing  problems  such  as  photosynthesis, 
ecosystem  dynamics,  ocean  health,  seawater  clarity,  underwater  imaging,  biogeochemical 
cycling,  carbon  budgets,  upper-ocean  thennodynamics,  and  climate  change  [Dickey  et 
al,  2011], 

The  simplest  optical  interactions  occur  at  the  boundary  of  the  atmosphere  and 
ocean,  and  are  governed  by  Snell’s  law  and  the  Fresnel  equations.  In  order  to  accurately 
model  the  very  complex  interactions  between  water  molecules  and  other  constituents, 
radiative-transfer  methods  and  precise  measurements  are  required  [Dickey  et  al,  201 1], 

Optical  remote  sensing  is  an  important  tool  utilized  for  monitoring  marine 
environments.  Since  light  is  readily  absorbed  by  water,  remote  sensing  is  most 
successful  in  places  with  very  shallow,  clear  water,  usually  up  to  depths  of  30  or  possibly 
40  meters.  Sea  water  contains  dissolved  and  particulate  matter,  with  varying 
concentrations  (both  spatially  and  temporally)  throughout  the  water  column  [ Mobley , 
1994;  Camacho,  2006].  Optical  properties  of  the  water  column  can  be  divided  into  two 
classes — inherent  optical  properties  (IOPs)  and  apparent  optical  properties  (AOPs)  [Smith 
and  Baker,  1981;  Camacho,  2006].  These  classes  are  further  discussed  in  sections  to 
follow. 


1,  Inherent  Optical  Properties  (IOPs) 

An  optical  property  is  inherent  if  it  depends  only  upon  the  medium  and  is 
independent  of  the  ambient  light  field  within  that  medium  [ Mobley ,  1994;  Camacho, 
2006].  As  light  enters  the  water  column,  it  interacts  with  particles  and  sediment,  which 
cause  the  incident  light  to  be  altered  by  scattering  or  absorption  [ Thomas  and  Stamnes, 
1999;  Camacho,  2006].  It  is  these  scattering  and  absorption  characteristics  that  define 
IOPs  of  water  [Camacho,  2006]. 
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2.  Apparent  Optical  Properties  (AOPs) 

An  optical  property  is  apparent  if  it  is  dependent  on  the  medium,  as  well  as  the 
directional  structure  of  the  ambient  light  field.  Similar  to  IOPs,  AOPs  are  also  dependent 
on  the  dissolved  particles  and  sediment  in  the  water  column.  Unlike  IOPs,  however, 
these  properties  cannot  be  measured  in  situ  because  they  depend  on  the  ambient  radiance 
[Mobley,  1994;  Camacho,  2006], 

3.  Sun  Glint  Correction  for  Spectral  Imagery 

Sun  glint  is  the  specular  reflection  of  light  directly  transmitted  from  the  upper  side 
of  the  air- water  interface  [Kay  et  ah,  2009].  It  typically  forms  bands  of  white  along  wave 
edges  on  the  windward  side  of  nearshore  environments  [Hed/ey  et  al.,  2005],  As  it  is  a 
serious  confounding  factor  for  remote  sensing  of  water  column  properties  and  benthos, 
researchers  have  spent  a  good  deal  of  time  creating  and  testing  techniques  to  estimate  and 
remove  the  glint  radiance  component  from  imagery.  Sun  glint,  which  is  a  function  of  sea 
surface  state,  sun  position,  and  viewing  angle,  occurs  in  imagery  when  the  water  surface 
orientation  is  such  that  the  sun  is  directly  reflected  toward  the  sensor.  The  component  of 
sensor-received  radiance  can  be  so  high  that  it  becomes  impossible  to  retrieve  any 
information  about  the  ocean  environment  [Kay  et  al.,  2009].  Examples  of  sun  glint 
appearing  in  imagery  are  shown  in  Figures  1 1  and  12. 
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NOAA  GOES  12 


091202  0845  UTC 


NASA  GSFC  GOES  Project 


Figure  1 1 .  An  example  of  sun  glint  on  the  ocean’s  surface  (image  captured  by  the 

Geostationary  Operational  Environmental  Satellite  (GOES)  on  2  December  2009) 
(Image  from  http://rsd.gsfc.nasa.gov/goes/text/hotstuff.html) 


Summer  Sunglint  in  the  Sargasso  Sea 
dawns  on  Cuba  and  Florida 
1200  GMT  on  22  June  2000 
GOES-West  visible  image 


GOES  project  NASA-GSFC 


NASA  Q 


Figure  12.  Another  example  demonstrating  the  effect  of  sun  glint  on  the  imager 
(image  captured  by  the  GOES  satellite  on  22  June  2000)  (Image  from 
http  ://rsd.  gsfc  .nasa.  gov/ goes/text/ho  tstuff.html) 
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Sun  glint  causes  problems  within  the  imagery  domain.  A  case  study  performed 
by  Goodman  et  al.  [2008]  reported  that  uncorrected  glint  in  high  resolution  imagery  led 
to  errors  as  large  as  30%  when  measuring  ocean  depth.  For  airborne  surveys,  optimal 
flight  paths  and  directions  can  be  chosen  based  on  the  time  of  day,  but  this  is  more 
difficult  in  the  case  of  satellite  imagery  [ Goodman  et  al.,  2008]. 

Glint  correction  methods  have  been  developed  to  improve  image  accuracy  for  two 
main  categories  of  water  types:  open  ocean  imaging,  and  higher  resolution  coastal  and 
aerial  applications.  The  aim  for  both  cases  is  the  same — to  estimate  the  glint  contribution 
to  the  radiance  reaching  the  sensor,  and  then  subtract  it  from  the  received  signal  [Kay  et 
al,  2009], 

The  set  of  methods  utilized  for  coastal  images  is  meant  to  be  used  with  pixel  sizes 
of  less  than  approximately  10  meters.  Data  from  the  NIR  band  are  used  as  an  indication 
of  the  amount  of  glint  in  the  received  signal.  The  spectrum  from  a  section  of  the  image 
containing  deep  water  is  then  used  to  establish  the  relationship  between  the  NIR  and  glint 
radiances  [Hochberg  et  al.,  2003;  Hedley  et  al,  2005;  Lyzenga  et  al,  2006;  Goodman  et 
al,  2008;  Kay  et  al,  2009].  The  assumption  of  a  linear  relationship  between  NIR 
brightness  and  the  amount  of  sun  glint  in  the  VIS  bands  holds  because  the  real  index  of 
refraction  is  approximately  equal  for  NIR  and  VIS  wavelengths  [Mobley,  1994;  Hedley  et 
al.,  2005]. 

When  no  sun  glint  is  present,  the  radiance  received  by  a  satellite-borne  sensor  is 
dominated  by  atmospheric  scattering.  Light  from  paths  through  the  air  constitute  over 
80%  of  received  radiance,  water-leaving  paths  make  up  approximately  15%,  and 
reflected  light  only  1  to  2%  [Sturm,  1981;  Kay  et  al.,  2009].  Sun  glint  can  increase  the 
reflected  radiance  by  a  factor  of  2  or  more,  and  the  worst  case  glint  can  saturate  the 
sensor,  making  it  impossible  to  retrieve  the  water-leaving  radiance  for  those  pixels  [Kay 
et  al.,  2009]. 
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D.  MULTI-ANGLE  RESEARCH 

Multi-angle  remote  sensing  capabilities  offer  a  number  of  advantages  with  respect 
to  a  single  shot  dataset.  Multi-angular  data  fusion  has  been  shown  to  allow: 

•  The  exploitation/investigation  of  BRDF, 

•  The  extraction  of  digital  height  maps  (DHMs), 

•  Atmospheric  parameter  retrieval, 

•  Classification  improvement,  etc.  [Pacifici  et  al.,  201 1] 

The  implications  of  adding  height  data  and  multi-angle  MS  reflectance,  both 
derived  from  the  multi-angle  sequence,  to  the  textural,  morphological,  and  spectral 
information  of  a  single  MS  image  has  been  investigated  by  Longbotham  et  al.  [2011].  It 
was  detennined  that  multi-angle  collections  significantly  increase  the  dimensionality  of 
the  data  available  for  a  single  target,  as  well  as  allow  for  differentiation  of  classes  not 
typically  well  identified  by  a  single  image  [ Longbotham ,  Chaapel  et  al,  201 1]. 

For  the  particular  Longbotham  et  al.  [2011]  study  referred  to  here,  the  analyzed 
data  sequence  was  collected  over  urban  areas  of  Atlanta,  GA  in  December  2009.  The  full 
dataset  contained  27  images,  with  satellite  elevation  (the  angle  between  the  horizon  and 
satellite  as  viewed  from  the  image  target)  as  low  as  25  degrees,  but  the  study  only  utilized 
those  images  with  a  moderate  (relative  to  the  satellite’s  capability)  off-nadir  observation 
angle  of  less  than  30  degrees.  This  left  the  group  with  13  images  between  57  degrees 
(forward-looking)  to  81.5  degrees  (most  nadir)  to  59  degrees  (backward-looking). 
Results  showed  a  27%  improvement  in  classification  accuracy  for  spatial  experiments, 
and  also  a  14%  improvement  in  classification  accuracy  for  spectral  experiments.  Figure 
13  shows  the  ground  observed  azimuth  and  elevation  of  the  WV-2  satellite  for  each 
image  observation,  as  well  as  the  azimuth  and  elevation  of  the  sun  during  the  multi-angle 
sequence  used  for  the  research  [Longbotham,  Chaapel  et  al.,  2011].  Further  sections  will 
discuss  why  WV-2  is  capable  of  acquiring  multi-angle  imagery. 
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Figure  13.  Ground  observed  azimuth  and  elevation  of  the  WV-2  satellite 

(black  markers)  for  each  image  observation  as  well  as  the  azimuth  and  elevation 
of  the  sun  (yellow  circle)  during  the  2009  Atlanta,  GA  multi-angle  sequence 
acquisition.  Azimuth  is  plotted  angularly  clockwise  (North:  0  degrees,  East:  90 
degrees,  South:  180  degrees,  West:  270  degrees),  and  elevation  is  plotted  radially 
from  the  center  (ground  nadir:  90  degrees,  ground  horizon:  0  degrees)  (From 
Longbotham,  Chaapel  et  al.  [2011]) 


There  is  significant  improvement  shown  over  the  baseline  classification  when 
using  a  multi-angle  WV-2  sequence.  In  terms  of  spatial  classification,  Longbotham  et  al. 
[2011]  successfully  demonstrated  the  ability  to  differentiate  between  classes  like  bridges 
and  man-made  structures,  which  are  generally  difficult  to  classify  because  they  are 
spectrally  similar  to  ground-level  classes  of  the  same  material  [ Longbotham  et  al,  2011]. 
Improvements  were  also  made  in  terms  of  spectral  classification — the  group  was  able  to 
distinguish  between  classes  that  may  prove  to  be  valuable  in  land-use  classification,  such 
as  moving  versus  parked  vehicles.  The  ability  to  differentiate  between  spectrally  similar 
classes  such  as  paved  parking  lots  and  paved  highways  were  also  made  possible  using  the 
multi-angle  sequence  [Longbotham,  Bleiler  et  al.,  2011]. 
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E. 


PREVIOUS  WORK  AT  THE  NAVAL  POSTGRADUATE  SCHOOL 


1.  Work  in  Satellite  Derived  Bathymetry 

a.  Depth  Analysis  of  Midway  Atoll  Using  QuickBird  Multi-Spectral 
Imaging  over  Variable  Substrates  (Camacho,  2006) 

Camacho  [2006]  used  QuickBird  MSI  over  Midway  Atoll,  Northwestern 
Hawaiian  Islands  in  an  attempt  to  analyze  depth  and  identify  variable  bottom-types  in 
shallow  water.  The  objective  of  this  study  was  to  use  an  MS  image  to  categorize  benthic 
substrates  based  on  spectral  characteristics  and  ground  truth  data  collected  in  situ.  With 
this  image,  a  ratio  of  reflectances  (the  “ratio  method”)  was  used  to  extract  depth 
separately  over  variable  substrates.  The  motivation  for  this  work  originated  from  a 
limitation  pointed  out  by  another  Naval  Postgraduate  School  thesis  student,  Clark  [2005], 
whose  results  demonstrated  that  the  ratio  method  proved  to  be  sensitive  to  bottom  type. 
The  method  produced  shallower  depths  over  bottom  types  with  low  albedo  and  deeper 
depths  over  bottom  types  with  high  albedo.  Clark  [2005]  also  noted  that  sun  glint  had  an 
effect  on  the  overall  results.  The  maximum  depth  that  could  be  obtained  using  the  ratio 
method  was  15  meters  [Clark,  2005]. 

Camacho  [2006]  utilized  methods  by  Lyzenga  [1978]  and  Stumpf  et  al. 
[2003]  to  extract  water  depth  and  bottom  type  infonnation,  and  also  used  the  ratio 
method  to  retrieve  accurate  depths  over  variable  bottom  types  and  low-albedo 
environments  (respectively)  [Lyzenga,  1978;  Stumpf  et  al.,  2003;  Camacho,  2006], 

b.  Depth  Derivation  from  the  WorldView-2  Satellite  Using 
Hyperspectral  Imagery  (Loomis,  2009) 

Loomis  [2009]  created  a  simulation  before  WV-2’s  launch  to  determine 
the  usefulness  of  the  Yellow  band  for  depth  derivation.  Loomis  [2009]  utilized  data  from 
the  Airborne  Visible/InfraRed  Imaging  Spectrometer  (AVIRIS)  sensor  over  Kaneohe 
Bay,  Hawaii,  and  then  processed  the  imagery  using  the  Stumpf  et  al.  [2003]  ratio  method 
to  detennine  bathymetry.  The  Green  &  Blue,  Yellow  &  Green,  and  Yellow  &  Blue  band 
combinations  were  compared  to  ground  truth  bathymetry  derived  from  a  digital  nautical 
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chart.  Results  indicated  that  use  of  the  Yellow  band  improved  the  accuracy  of  derived 
depths,  especially  in  shallow  water  [Loomis,  2009]. 

c.  Contributions  to  Remote  Sensing  of  Shallow  Water  Depth  with 
the  WorldView-2  Yellow  Band  (Madden,  2011) 

Following  Loomis’  [2009]  work,  Madden  [2011]  analyzed  bathymetry  in 
Tampa  Bay,  Florida’s  shallow  water  using  WV-2  imagery.  The  Yellow  band  was 
combined  separately  with  the  Blue,  Green,  and  Red  bands,  and  then  compared  to  the 
more  traditional  Blue  &  Green  and  Green  &  Red  combinations.  Madden ’s  [2011]  results 
showed  that  the  addition  of  the  Yellow  band  provided  more  information  about 
bathymetry,  but  less  sensitivity  to  bottom  type  in  two  of  three  transect  lines  used 
[Madden,  2011], 

2.  Work  in  Multi-Observation  Imagery 

a.  Coastal  Bathymetry  Using  8-Color  Multispectral  Satellite 
Observation  of  Wave  Motion  (McCarthy,  2010) 

McCarthy  [2010]  measured  coastal  bathymetry  near  Camp  Pendleton  in 
California  by  using  wave  motion  as  observed  by  WV-2.  After  principal  component 
transforms  were  performed,  successive  change  detection  images  contained  both  spatial 
and  temporal  information.  Wave  celerity  could  be  detennined  and  depth  inversion  was 
performed.  McCarthy  [2010]  measured  the  wavelength  of  a  wave  above  a  known  depth, 
and  then  used  the  wave  period  method  to  detennine  depth  for  other  waves  in  the 
propagation  direction  of  that  particular  wave.  McCarthy  [2010]  detennined  that  the 
spatial  resolution  for  this  method  was  higher  and  possibly  more  accurate  than  the 
reference  bathymetry  used,  especially  in  the  surf  zone  [McCarthy,  2010], 

b.  Coastal  Bathymetry  Using  Satellite  Observation  in  Support  of 
Intelligence  Preparations  of  the  Environment  (Myrick,  2011) 

Myrick  [2011]  followed  McCarthy’s  [2010]  work,  and  calculated  water 
depths  using  linear  dispersion  relationships  for  surface  gravity  waves.  Myrick  [201 1]  was 
able  to  establish  depth  profiles  out  to  1  kilometer  offshore,  as  well  as  derive  depths  up  to 
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15  meters.  Comparisons  with  United  States  Geological  Survey  (USGS)  bathymetric 
acquisitions  from  2009  show  agreement  within  5%  in  the  surf  zone  and  1%  outside  of  the 
surf  zone  [ Myrick ,  2011]. 

c.  Automating  Nearshore  Bathymetry  Extraction  from  Wave 
Motion  in  Satellite  Optical  Imagery  (Mancini,  2012) 

Mancini  [2012]  extracted  nearshore  depths  for  Waimanalo  Beach,  Hawaii 
from  WV-2  optical  imagery  by  means  of  automated  wave  kinematics  bathymetry  (WKB). 
Two  sets  of  three  sequential  images,  acquired  at  approximately  10  second  intervals,  were 
used.  Depths  from  close  to  shore  out  to  about  20  meters  deep  were  generated. 
Comparisons  to  the  Scanning  Hydrographic  Operational  Airborne  Lidar  Survey 
(SHOALS)  LiDAR  bathymetry  values  showed  WKB  depths  were  accurate  to  about  0.5 
meters,  with  R  values  of  90%,  and  were  frequently  in  the  range  of  10%  to  20%  relative 
error  for  depths  ranging  from  2  to  16  meters  [ Mancini ,  2012], 

d.  High  Spatial  Resolution  Bidirectional  Reflectance  Retrieval 
Using  Satellite  Data  (McConnon,  2010) 

McConnon  [2010]  analyzed  15  WV-2  images  acquired  over  Duck,  North 
Carolina  and  10  images  of  Pendleton,  California  to  extract  BRDF.  Separation  angles 
between  the  inbound,  solar  zenith  angle,  and  the  outbound  reflectance  angle  were 
calculated.  BRDF  shifts  across  wavelengths  and  regions  of  interest  (ROIs)  were  plotted 
against  the  previously  mentioned  values  and  then  examined  [McConnon,  2010]. 
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III.  DATA  AND  METHODS 


A.  DATA 

1.  Satellite  Sensor,  Study  Site,  and  Imagery  Dataset 
a.  WorldView-2  (WV-2) 

WorldView-2,  launched  on  8  October  2009  from  Vandenberg  Air  Force 
Base,  is  DigitalGlobe’s  third  operational  satellite  after  QuickBird  and  WorldView-1.  The 
sensor  is  in  a  nearly  circular,  sun-synchronous  orbit,  and  flies  at  an  altitude  of 
approximately  770  kilometers  (Table  1).  At  nadir,  the  best  possible  spatial  resolution  is 
0.46  meters  (panchromatic)  and  1.84  meters  (MS).  What  differentiates  WV-2  from 
DigitalGlobe’s  previous  sensors  is  its  9  spectral  bands — one  panchromatic  (PAN), 
ranging  from  approximately  450  to  800  nm  (centered  at  632  nm),  and  8  MS,  ranging 
from  approximately  400  to  1050  nm  (Figure  14,  Table  2).  These  MS  bands  include: 
Coastal  (centered  at  427  nm),  Blue  (centered  at  478  nm),  Green  (centered  at  546  nm), 
Yellow  (centered  at  608  nm),  Red  (centered  at  659  nm),  Red  Edge  (centered  at  724  nm), 
NIR-1  (centered  at  831  nm),  and  NIR-2  (centered  at  908  nm)  (available  at 
http://www.digitalglobe.com/about-us/content-collection#satellites&worldview-2). 
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Table  1.  WV-2  Design  and  Specifications  (After  infonnation  available  at 
http://www.digitalglobe.com/about-us/content- 


collection#satellites&worldview-2) 

Launch  Information 

Date:  October  8,  2009 

Launch  Site:  Vandenberg  AFB,  CA,  USA 

Orbit 

Altitude:  770  km 

Type:  Sun  synchronous 

Period:  100  minutes 

Sensor  Resolution 

Panchromatic:  0.46  m  GSD  at  nadir,  0.52  m  GSD  at  20 

degrees  off-nadir 

Multispectral:  1.85  m  GSD  at  nadir,  2.07  m  GSD  at  20 

degrees  off-nadir 

Swath  Width 

16.4  km  at  nadir 

Attitude  Determination 

3 -axis  stabilized 

and  Control 

Actuators:  Control  Moment  Gyros  (CMGs) 

Sensors:  Star  trackers,  solid  state  IRU,  GPS 

Pointing  Accuracy  and 

Accuracy:  <500  m  at  image  start  and  stop 

Knowledge 

Knowledge:  Supports  geolocation  accuracy  below 

Retargeting  Agility 

Time  to  Slew  200  km:  10  seconds 

Revisit  Frequency  (at 

1 . 1  days  at  1  m  GSD  or  less 

40  degrees  N  Latitude) 

3.7  days  at  20  degrees  off-nadir  or  less  (0.52  m  GSD) 

Geolocation  Accuracy 

Demonstrated  <3.5  m  CE90  without  ground  control 

Capacity 

1  million  sq  km  per  day 
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WV02  Relative  Spectral  Radiance  Response 


Figure  14.  WV-2’s  relative  spectral  radiance  response 
(From  Updike  and  Comp  [2010]) 

Table  2.  WV-2  Sensor  Bands  (After  infonnation  available  at 
http://www.digitalglobe.com/about-us/content- 
collection#satellites&worldview-2) 


Band 

Band  Width  (nm)  Center  Wavelength  (nm) 

Panchromatic  450  -  800 

632 

Coastal 

400-450 

427 

Blue 

450-510 

478 

Green 

510-580 

546 

Yellow 

585  -625 

608 

Red 

630  -  690 

659 

Red  Edge 

705  -  745 

724 

NIR-1 

770  -  895 

831 

NIR-2 

860  -  1040 

908 
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The  WV-2  system’s  8  MS  bands  are  arranged  in  two  arrays  of  4  MS  bands 
each  (MSI  and  MS2).  MSI  includes  Blue,  Green,  Red,  and  NIR-1.  MS2  is  comprised  of 
Coastal,  Yellow,  Red  Edge,  and  NIR-2.  Imaging  options  are:  PAN  only,  PAN  +  MSI, 
and  PAN  +  8  MS  (MSI  and  MS2).  There  is  a  0.3  second  delay  between  MSI  and  MS2 
acquisitions,  according  to  G.  Miecznik  (unpublished  data,  2012). 

When  combined,  these  bands  are  designed  to  improve  the  segmentation 
and  classification  of  land  and  aquatic  features  beyond  any  other  multispectral  satellite 
imager.  After  WV-2’s  launch,  it  was  speculated  that  the  increased  agility  and  addition  of 
the  Coastal  band  would  improve  remote  bathymetric  measurements  (mainly  due  to  the 
Coastal  band’s  wavelength  value,  making  it  least  absorptive  by  water).  Analysts 
expected  to  be  able  to  calculate  depths  up  to  20  meters,  and  possibly  even  30  meters 
using  the  Coastal,  Blue,  and  Green  bands.  Once  scientists  were  able  to  utilize  WV-2 
imagery,  the  Coastal  band  has  proven  to  be  useful  for  the  retrieval  of  water  depth,  true- 
color  correction  for  human  vision  representation,  chlorophyll  absorption,  and 
atmospheric  scattering  correction  [Pacifici  and  Navular,  2011;  Marchisio  et  al.,  2011]. 
Also,  with  an  average  revisit  time  of  1.1  days  and  the  ability  to  utilize 
off-nadir  pointing,  WV-2  was  expected  to  better  monitor  ocean  environments,  which 
are  highly  dynamic  and  constantly  changing  (information  available  at 
http://www.satimagingcorp.com/media/pdf/Bathymetry_Datasheet.pdf). 

Bi-directional  scanning  is  supported  by  WV-2.  The  onboard  camera  has  a 
standard  maximum  look  angle  of  approximately  40  degrees  off-nadir  and  can  slew  across 
300  kilometers  of  the  Earth’s  surface  in  9  seconds.  This  unique  agility  enables  WV-2  to 
collect  dense  image  sequences  of  a  single  target  from  several  observation  angles  as  it 
overflies  an  area  [Longbotham  et  al.,  201 1]. 

b.  Kailua  Bay,  Oahu,  Hawaii 

This  study  focused  on  one  location  on  the  windward  side  of  the  Oahu 
coastline  of  Hawaii  -  Kailua  Bay  (Figure  15).  The  bay’s  approximate  latitude/longitude 
is:  21  degrees  24’  29”  N,  157  degrees  44’  09”  W.  This  particular  beach  was  chosen 
because  it  was  the  least  cloudy  out  of  the  39  WV-2  images  in  the  acquired  dataset. 
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Figure  15.  The  Hawaiian  Islands  (left),  focusing  on  Kailua  Bay,  Oahu  (right)  (From 

Google  Earth) 


Kailua  Bay  is  a  carbonate  reef-dominated  embayment.  There  are  two 
categories  of  benthic  substrate  found  here:  areas  of  carbonate  sand  and  fossil  reef 
hardgrounds,  and  reef  habitats  of  coral  and  algae  species.  There  is  a  sand-floored  channel 
at  the  center  of  the  bay,  which  cuts  across  the  reef  and  connects  the  seaward  and 
nearshore  sand  fields.  Algae  and  corals  grow  on  the  plains  [ Isoun  et  al.,  2003]. 


Areas  with  sand  and  fossil  reef  appear  light-colored  and  are  highly 
reflective  in  the  WV-2  imagery.  The  coral  and  algae  communities  look  dark  and  have 
low  reflectance  [Isoun  et  al.,  2003].  Aerial  images  are  shown  in  Figure  16. 


Figure  16.  Oblique  aerial  imagery  acquired  over  Kailua  Bay  on  10  November  2003 
on  a  Nikon  Coolpix  5700  digital  camera  (settings  and  focus  all  automated) 
from  a  Cessna  plane  flying  at  about  2,000  feet  (From 
http://www.soest.hawaii.edu/coasts/erosion/oahu/oblique.php) 
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c.  Imagery  Dataset 

A  collection  of  39  multispectral  images  of  the  windward  side  of  Oahu, 
Hawaii  was  acquired  by  the  WorldView-2  satellite  on  30  July  2011  between  21:22:49Z 
and  21:28:54Z.  Images  were  acquired  at  approximately  10  second  intervals,  covering 
just  over  6  minutes  (Figure  17).  Section  A  of  the  Appendix  contains  a  sample  WV-2 
metadata  file. 


Figure  17.  Five  examples  out  of  the  39  WV-2  image  acquisitions 

(labels  are  actual  Image  IDs) 


2.  Bathymetry  Data 

Bathymetry  data  of  Kailua  Bay  were  acquired  from  the  University  of  Hawaii  at 
Manoa,  School  of  Ocean  &  Earth  Science  &  Technology,  Department  of  Geology  and 
Geophysics  (data  are  available  at  http://www.soest.hawaii.edu/coasts/ 
data/oahu/shoals. html  and  information  about  this  dataset  is  available  at 
http://www.soest.hawaii.edU/coasts/data/readme.html#shoals).  According  to  the 
University  of  Hawaii  website,  these  bathymetry  data  were  collected  from  the  SHOALS 
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website,  as  part  of  a  survey  conducted  in  2000  (information  is  available  at 
http://shoals.sam.usace.army.mil/hawaii/pages/Oahu.htm).  After  investigating  the  origins 
of  the  dataset,  however,  and  determining  that  these  data  were  not  collected  by  SHOALS 
in  2000,  it  is  now  believed  that  the  data  were  collected  by  USGS  circa  2002  to  2005. 
This  information  was  provided  by  C.  Fletcher  (unpublished  data,  2012).  Prior  to  the 
completion  of  this  particular  project,  there  was  no  resolution  concerning  the  actual  details 
about  the  dataset. 

3.  Software 

a.  Environment  for  Visualizing  Images  +  Interactive  Data 
Language  4.8  (ENVI  4.8  +  IDL) 

Environment  for  Visualizing  Images  (ENVI)  4.8  was  the  main  software 
program  used  for  data  analysis  and  manipulation.  ENVI  is  an  image  processing  software 
system  designed  for  multi-  and  hyperspectral  data  analysis  and  information  extraction. 
This  software  is  written  in  the  Interactive  Data  Language  (IDL),  a  programming  language 
that  provides  integrated  image  processing  and  display  capabilities  [Research  Systems, 
2004;  Camacho,  2006]. 

ENVI  was  used  to  process  the  WV-2  imagery.  All  image  mosaicking, 
spatial  subsetting  and  chipping,  radiance  calibration,  and  registration  was  performed  in 
ENVI.  A  number  of  IDL  programs  were  written  for  each  combination  of  two  WV-2 
bands  tested,  and  each  incorporated  a  section  to  create  masks  for  land,  clouds,  whitecaps, 
and  glint,  as  well  as  apply  a  band  ratio  method  for  depth  determination  (see  Section  C  of 
the  Appendix). 

B.  METHODS 

A  number  of  pre-processing  steps  were  perfonned  before  bathymetric  derivation 
could  occur.  First,  images  were  mosaicked  and  analyzed  for  overall  quality  and  cloud 
cover.  The  next  steps  included  map  coordinate  conversion,  radiance  calibration, 
land/cloud  mask  creation,  sun  glint  removal,  and  the  application  of  a  ratio  method  for 
bathymetry  derivation. 
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1.  Bathymetry  Derivation 

a.  Assess  Data  Coverage,  Angles,  and  Quality 

Before  any  data  preparation  steps  occurred,  images  were  assessed.  It  was 
important  to  ensure  that  imagery  were  collected  over  the  correct  site,  and  that  datasets  did 
not  include  any  bad  data.  Scenes  were  ordered  by  mean  satellite  elevation  angle  (from 
most  forward-looking,  to  most  nadir,  to  most  backward-looking).  Table  3  and  Section  B 
of  the  Appendix  list  all  39  images  acquired  over  Oahu. 
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Table  3. 


firstLineTime 

(Z) 

meanSatEl 

meanSatAz 

21:22:49 

24.9 

18.8 

21:22:58 

26.5 

19.2 

21:23:07 

28.1 

19.6 

21:23:16 

29.8 

20.1 

21:23:25 

31.7 

20.6 

21:23:34 

33.7 

21.2 

21:23:53 

38.2 

22.6 

21:24:02 

40.8 

23.5 

21:24:12 

43.7 

24.6 

21:24:21 

46.7 

26.0 

21:24:31 

50.1 

27.6 

21:24:41 

53.6 

29.7 

21:24:51 

57.6 

32.4 

21:25:01 

61.6 

36.1 

21:25:12 

66.0 

41.5 

21:25:22 

70.2 

49.2 

21:25:32 

74.0 

61.2 

21:25:42 

76.9 

80.1 

21:25:52 

77.8 

105.5 

21:26:02 

76.3 

129.8 

21:26:12 

73.1 

146.9 

21:26:23 

69.0 

157.9 

21:26:33 

64.8 

164.7 

firstLineTime 


meanSatEl 

meanSatAz 

21:26:43 

60.6 

169.4 

21:26:53 

56.6 

172.8 

21:27:02 

52.9 

175.3 

21:27:12 

49.6 

177.2 

21:27:21 

46.5 

178.7 

21:27:30 

43.7 

179.9 

21:27:39 

41.1 

180.9 

21:27:47 

38.7 

181.7 

21:27:56 

36.5 

182.5 

21:28:04 

34.4 

183.1 

21:28:13 

32.5 

183.7 

21:28:21 

30.7 

184.2 

21:28:29 

29.0 

184.6 

21:28:38 

27.4 

185.0 

21:28:46 

25.9 

185.4 

21:28:54 

24.5 

185.7 

b.  Mosaic  Data 

Data  files  were  delivered  in  such  a  way  that  scenes  needed  to  be 
mosaicked.  This  step  was  performed  to  be  able  to  better  view  the  entire  coastline. 
Spatial  subsetting  could  then  occur  without  worrying  about  areas  that  may  not  have  full 
coverage  because  the  rest  of  the  scene  was  in  a  different  file.  Figure  18  shows  example 
images  after  mosaicking. 


Most  Forward-Looking  Most  Nadir  Most  Backward-Looking 

(24.9  degrees)  (77.8  degrees)  (24.9  degrees) 


Figure  18.  The  most  forward-looking,  nadir,  and  backward-looking  images 

after  data  coverage  assessment  and  mosaicking 

c.  Convert  Map  Projection 

Level  IB  (Basic)  images  were  used,  meaning  that  data  were  not  projected 
to  a  plane  using  a  map  projection  or  datum  and  needed  to  be  converted  to  Universal 
Transverse  Mercator  (UTM),  World  Geodetic  System  1984  (WGS-84),  Zone  4N 
(information  is  available  at  http://www.digitalglobe.com/downloads/ 
DigitalGlobe_Core_Imagery_Products_Guide.pdf).  Pixels  were  also  made  square. 
Figure  19  shows  how  the  map  projection  conversion  affected  the  look  of  the  imagery. 
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Figure  19.  Map  projection  conversion  from  Geographic  Latitude/Longitude  (left) 

to  UTM,  WGS-84,  Zone  4N  (right) 


d.  Radiance  Calibration 

The  WV-2  spectral  radiance  response  is  defined  as  the  ratio  of  the  number 
of  photo-electrons  measured  by  the  system,  to  the  spectral  radiance  [W-rn  -sr  -pm  ]  at  a 
certain  wavelength  present  at  the  entrance  to  the  telescope  aperture.  The  spectral 
radiance  response  for  each  band  is  normalized  by  dividing  by  the  maximum  response 
value  for  that  band  to  arrive  at  a  relative  spectral  radiance  response  [ Updike  and  Comp, 
2010]. 

Relative  radiometric  calibration  and  correction  are  necessary.  This  is 
because  a  uniform  scene  does  not  create  a  uniform  image  when  it  comes  to  raw  digital 
numbers  (DNs).  This  type  of  correction  minimizes  image  artifacts,  such  as  vertical 
streaks  or  bands  due  to  differences  in  gain  or  offset,  and  is  performed  on  raw  data  from 
all  detectors  in  all  bands  during  the  early  stages  of  WV-2  product  generation.  The 
products  are  linearly  scaled  to  absolute  spectral  radiance  [ Updike  and  Comp,  2010]. 

In  the  case  of  large  mosaics,  radiometric  balancing  will  help  match  the 
brightness  of  the  other  scenes  used  in  the  mosaic  [ Updike  and  Comp,  2010].  As  glint 
removal  is  performed  on  images  after  the  radiance  calibration  step,  all  images  were  left  in 
radiance  for  this  research. 
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e.  Subset  Data 

Images  were  subset  to  focus  on  Kailua  Bay  (Figure  20).  This  was  done  to 
reduce  file  size,  emphasize  a  spot  with  interesting  and  variable  bathymetry,  as  well  as  to 
better  concentrate  on  a  specific  location  that  had  less  cloud  cover. 


Figure  20.  The  full  scene  (left)  was  chipped  to  focus  on  Kailua  Bay  (right) 

f  Registration  of  Off-Nadir  Images  to  the  Most  Nadir  Image 

All  non-nadir  images  were  registered  to  the  most  nadir  image.  Twenty  tie 
points  were  interactively  selected  for  the  same  locations  in  all  of  the  non-nadir  images. 
The  data  were  then  warped  using  a  first  order  polynomial.  Maximum  pixel  error  for  the 
registrations  was  1.521  pixels  with  an  average  root  mean  square  (RMS)  error  of  about 
0.895  pixels  for  the  entire  dataset.  Images  were  then  re-chipped  so  that  every  scene 
covered  the  same  geographic  coordinates.  All  chips  were  reduced  to  995  samples  and 
999  lines.  An  example  is  shown  in  Figure  2 1 . 
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1010  2010  2100  3100  4100 


Figure  2 1 .  After  registration,  each  image  (top)  was  chipped  to  the 

995  x  999  pixel  scene  shown  (bottom) 


g.  Land,  Glint,  Cloud,  and  Whitecap  Masks 

Application  of  spatial  masks  for  land,  glint,  cloud,  and  whitecap  areas 
allows  analysis  to  be  limited  to  only  the  areas  and  materials  of  interest.  Masks  were 
determined  by  comparing  scatter  plots  using  the  Blue  and  NIR-1  bands.  Water  regions 
are  those  in  the  lower  left  portion  of  the  scatter  plot,  as  can  be  seen  in  Figure  22. 


Figure  22.  The  original  image  (left)  has  a  corresponding  scatter  plot  (middle); 

specific  sections  of  the  plot  are  highlighted  to  show  how  they  correspond  to  the 
image  (right) — red  is  mostly  land,  green  is  mostly  whitecaps  and  sun  glint, 

and  blue  is  water 
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The  scatter  plot  classified  land,  glint,  clouds,  and  whitecaps  as  any  pixels 
with  values  greater  than  the  user-defined  points  chosen  within  the  plot.  Every  other  pixel 
would,  therefore,  be  considered  not  land,  glint,  clouds,  or  whitecaps  and  would  not  be 
masked  (Figure  23).  Figure  24  illustrates  how  Figure  23 ’s  user-defined  region  affects  the 
imagery. 


2D  Scatter  Plot  Used  for  Glint  Removal 


6  8  10  12  14  16 

Blue  Band 


Figure  23.  2D  scatter  plot  of  the  Blue  vs.  NIR-1  bands;  the  user-defined  outline 

region  includes  pixels  containing  water,  and  everything  else  can  generally  be 
classified  as  land,  clouds,  glint,  whitecaps,  etc.  (e.g.,  not  water) 
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Figure  24.  Land,  glint,  clouds,  and  whitecaps  have  been  masked  in  this  image 

h.  Glint  Removal 

As  discussed  previously,  if  the  ocean  surface  is  perfectly  flat,  the 
reflection  of  the  sun  appears  as  a  bright,  small  portion  of  the  sea  surface.  If  the  ocean 
surface  is  not  flat  due  to  stronger  winds,  however,  parts  of  the  water  surface  further  from 
the  center  of  the  sun  glint  pattern  will  be  at  the  required  orientation  to  reflect  sun  light  to 
the  viewer.  Therefore,  the  sun’s  reflection  will  cover  a  larger  part  of  the  ocean  surface, 
and  will  be  made  up  of  many  tiny  highlights  that  each  reflect  from  a  particular  point  on 
the  surface  from  the  sun.  The  brightness  of  each  pixel  equals  the  total  brightness  of  all 
highlights  in  that  pixel  and  is,  therefore,  proportional  to  the  fraction  of  the  sea  surface  at 
the  right  slope  [ Kay  et  al.,  2009]. 

If  the  sun  and  sensor  are  treated  as  a  point  source  and  detector  for  a  given 
viewing  geometry,  there  is  only  one  facet  slope  and  orientation  that  is  consistent  with 
specular  reflection.  The  sun,  in  reality,  has  an  angular  diameter  of  0.53  degrees,  so  there 
is  a  range  of  possible  slopes  that  can  reflect  light  from  some  part  of  the  sun’s  disc  into  the 
sensor  [Kay  et  al.,  2009]. 

Steps  to  remove  glint  were  applied  to  the  two  WV-2  bands  that  were  the 
focus  of  each  version  of  code  (Section  C  of  the  Appendix).  The  method  used  was  based 
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on  the  method  revised  by  Hedley  et  al.  [2005]  after  Hochberg  et  al.  [2003],  Work  by 
Hedley  et  al.  [2005]  establishes  the  linear  relationship  between  NIR  and  VIS  bands  using 
a  linear  regression  based  on  a  sample  of  the  image  pixels.  Over  areas  with  underlying 
spectral  brightness,  such  as  deep  water,  one  or  more  regions  with  a  range  of  sun  glint  are 
selected.  For  each  VIS  band,  all  selected  pixels  are  included  in  a  linear  regression  of  NIR 
brightness  (x-axis)  against  the  VIS  band  brightness  (y-axis).  If  the  slope  of  this  line  for 
band  i  is  bj,  then  all  the  pixels  in  the  image  can  be  deglinted  in  band  i  by  applying  the 
following  equation: 

R  )  =  Ri  -  bj(R\\R  -  MinNiR),  (2) 

which  means:  reduce  the  pixel  value  in  band  i  ( R, )  by  the  product  of  regression  slope  (b,) 
and  the  difference  between  the  pixel  NIR  value  (Rnir)  and  the  ambient  NIR  level 
(MinNiR).  R  \  is  the  sun  glint  corrected  pixel  brightness  in  band  i.  MinNiR  represents  the 
NIR  brightness  of  a  pixel  with  zero  sun  glint,  and  can  be  estimated  by  the  minimum  NIR 
found  in  the  regression  sample  or  as  the  minimum  NIR  value  found  in  the  entire  image 
[Hedley  et  al.,  2005],  Refer  to  Figure  25  for  a  graphical  interpretation  of  the  Hedley  et  al. 
[2005]  method. 


Figure  25.  Graphical  interpretation  of  the  Hedley  et  al.  [2005]  de-glinting  method 

(From  Hedley  et  al.  [2005]) 
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The  analysis  incorporates  the  slope  of  the  regression  line.  A  simplified 
example  of  the  modified  equation  used  in  the  code  (Section  C  of  the  Appendix)  is  as 
follows: 

Green  =  Green  -  Slope  of  Green  Regression  Line  *  (NIR-1  -  MinNIR-1)  (3) 

Yellow  =  Yellow  -  Slope  of  Yellow  Regression  Line  *  (NIR-2  -  MinNIR-2)  (4) 

There  are,  however,  a  number  of  differences  between  the  Hedley  et  al.  [2005]  method 
and  the  one  used  for  this  research.  This  code  utilized  only  two  of  the  six  possible  WV-2 
VIS  bands  at  a  time.  It  also  focused  on  the  entire,  global  scene,  rather  than  a  small,  local 
portion  (only  deep  water,  for  example).  It  also  incorporated  masking,  which  was  not 
used  in  the  Hedley  et  al.  [2005]  research  [Hedley  et  al.,  2005]  (refer  to  Section  C  of  the 
Appendix). 


i.  Band  Ratio  to  Determine  Relative  Bathymetry 

A  two-step  process  was  used  to  derive  bathymetry.  Relative  bathymetry 
was  determined  by  performing  a  band  ratio  method,  and  then  derived  bathymetry  values 
were  obtained  by  regressing  relative  bathymetry  values  against  verified  depth  data  (to  be 
further  discussed  in  the  upcoming  sections).  Camacho  [2006]  calculated  relative 
bathymetry  using  the  natural  log  transformed  reflectance  values  of  the  deglinted 
reflectance  image  [ Camacho ,  2006]: 

ln(  1 000 *  b  1  )/ln(  1 000 *b2)  (5) 

The  relative  bathymetric  values  in  this  research  were  extracted  using  the  following 
expression,  with  a  slight  adjustment  to  the  constant  value  used  by  Camacho  [2006]: 

ln(  1 00  *  b  1  )/ln(  1 00  *b2)  (6) 

where  bl  is  the  first  band  used,  and  b2  is  the  second.  This  is  a  modification  of  the 
equation  used  by  Stumpf  et  al.  [2003]  [Stumpf  et  al.,  2003]. 
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Comparison  of  Water  Depths  to  LiDAR  Bathymetry 


The  bathymetry  data  of  Kailua  Bay  were  acquired  in  shapefile  format. 
They  were  then  converted  to  a  raster  in  the  proper  map  projection,  geographically  linked 
to  the  image  chip,  and  then  clipped  to  the  995  by  999  pixels-size.  These  data  were  then 
used  as  the  “true”  data. 

k.  Derived  Bathymetry 

Derived  depth  values  were  calculated  by  regressing  the  relative  depth 
values  with  the  actual  depth  values  collected  by  the  “true”  bathymetry.  The  chi-squared 
and  correlation  values  were  additional  IDL  code  outputs  that  helped  detennine  the 
“goodness  of  fit.”  A  lower  chi-squared  value  is  indicative  of  a  better  fit.  Correlation 
values  are  interpreted  as  percentages,  so  the  highest  value  represents  the  best  fit. 

2.  Analyze  Effects  of  Collection  Geometry  on  Water  Depth  Derivation 

The  effects  of  collection  geometry  on  water  depth  derivation  were  analyzed  by 
running  all  39  images  through  the  processing  and  analysis  approach  described  above. 
Chi-squared  and  correlation  values  were  recorded,  and  images  of  derived  depth  versus 
actual  “true”  depth  plots  were  compared.  Results  will  be  discussed  further  in  the 
following  section. 
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IV.  RESULTS 


A  total  of  39  bathymetric  maps  were  generated  from  the  WV-2  data  using  the 
previously  described  band  ratio  method.  The  accuracies  of  these  maps  were  then 
evaluated  using  “true”  bathymetry  data  from  the  University  of  Hawaii  at  Manoa.  The 
final  WV-2  derived  bathymetry  was  produced  by  regressing  the  relative  bathymetry 
values  against  the  actual  “true”  bathymetric  measurement.  The  chi-squared  and 
correlation  values  were  then  compared  to  analyze  effects  of  multi-angle  acquisition  on 
depth  derivation. 

A.  BATHYMETRY  FROM  ENTIRE  IMAGE 

A  total  of  39  bathymetric  maps  from  a  WV-2  multi-angle  sequence  over  Kailua 
Bay,  Oahu,  Hawaii  were  generated.  The  following  band  combinations  were  analyzed: 

•  Coastal  &  Blue, 

•  Coastal  &  Green, 

•  Blue  &  Green,  and 

•  Green  &  Yellow. 

Each  band  combination  was  analyzed  using  the  band  ratio  method  to  obtain  relative 
bathymetry.  These  values  were  then  regressed  against  the  “true”  bathymetry  data  shown 
in  Figure  26. 


43 


Figure  26.  “True”  bathymetry  data;  depth  is  represented  by  the  Rainbow  scale 

(shallow  water  is  red  and  deep  water  is  black) 

The  results  of  the  regression  were  the  derived  depth  values.  These  were  scaled  to 
meters  and  plotted.  Derived  depth  plots  from  Image  IDs  1010  (most  forward-looking), 
2010,  2100  (most  nadir),  3100,  and  4100  (most  backward-looking)  are  shown  in  Figures 
27  through  31  for  each  band  combination.  All  images  have  the  same  scale,  a  range  of 
depths  from  0  (red)  to  -18  (black)  meters. 
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1010 


Coastal  &  Blue  Coastal  &  Green  Blue  &  Green  Green  &  Yellow 

Figure  27.  Derived  depths  for  each  band  combination  for  Image  1010 

(most  forward-looking) 


2010 


Coastal  &  Blue  Coastal  &  Green  Blue  &  Green  Green  &  Yellow 

Figure  28.  Derived  depths  for  each  band  combination  for  Image  2010 


2100 


Coastal  &  Blue  Coastal  &  Green  Blue  &  Green  Green  &  Yellow 

Figure  29.  Derived  depths  for  each  band  combination  for  Image  2100 

(most  nadir) 
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3100 


Coastal  &  Blue  Coastal  &  Green  Blue  &  Green  Green  &  Yellow 


Figure  30.  Derived  depths  for  each  band  combination  for  Image  3100 


4100 


Coastal  &  Blue  Coastal  &  Green  Blue  &  Green  Green  &  Yellow 


Figure  3 1 .  Derived  depths  for  each  band  combination  for  Image  4100 

(most  backward-looking) 

It  can  be  seen  that  none  of  these  derived  depth  maps  come  close  to  a  perfect 
correlation  to  the  “true”  bathymetry  data.  This  most  likely  has  to  do  with  a  number  of 
factors,  including  high  off-nadir  acquisition  angle  and  also  cloud  cover  (seen  as  white, 
masked  data  in  Figures  30  and  31).  Due  to  the  fact  that  wave  patterns  can  be  seen  in  the 
derived  depth  images,  it  is  also  speculated  that  a  better  glint  removal  method  may  be 
required  in  order  to  derive  more  accurate  depths. 

Examination  of  these  images  shows  that  the  Coastal  &  Green  and  Blue  &  Green 
WV-2  band  combinations  tend  to  perform  better.  Also,  derived  depth  more  closely 
matches  the  “true”  depth  as  the  sensor  is  acquiring  images  closer  to  nadir. 
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B.  VARIATION  IN  ACQUISITION  ANGLES 

In  order  to  better  understand  the  effect  of  acquisition  angle  on  depth 
determination,  plots  of  derived  depth  versus  “true”  depth  for  Image  IDs  1010  (most 
forward-looking),  2010,  2100  (most  nadir),  3100,  and  4100  (most  backward-looking) 
were  created  and  compared.  The  chi-squared  and  correlation  values  of  each  were 
detennined  and  are  displayed  for  Image  ID  2100  (most  nadir)  in  Figures  32  through  35 
(the  rest  can  be  located  in  Section  D  of  the  Appendix). 


2100:  Coastal  Sc  Blue  Derived  vs.  True  Bathymetry 


-8  -7  -6  -5  -4  -3  -2 

Derived  Bathymetry  (m) 


Figure  32.  Coastal  &  Blue:  Derived  depth  vs.  “true”  depth  for  Image  2100 

(most  nadir) 
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2100:  Blue  Sc  Green  Derived  vs.  True  Bathymetry 


Derived  Bathymetry  (m) 


Figure  34.  Blue  &  Green:  Derived  depth  vs.  “true”  depth  for  Image  2100 

(most  nadir) 
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2100:  Green  &  Yellow  Derived  vs.  True  Bathymetry 


Derived  Bathymetry  (m) 


Figure  35.  Green  &  Yellow:  Derived  depth  vs.  “true”  depth  for  Image  2100 

(most  nadir) 

Chi-squared  and  correlation  values  were  calculated  and  plotted  for  each  Image  ID 
(view  angle)  for  all  four  band  combinations.  These  values  represent  the  correlation  of  the 
depth  derived  from  the  band  combination  ratios  versus  the  “true”  depth  (example  values 
are  shown  in  Section  E  of  the  Appendix).  Results  are  depicted  in  Figures  36  through  39. 
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Coastal  &  Blue  Chi-Squared 


Mean  Satellite  Elevation  Angle 

Figure  36.  Coastal  &  Blue:  Chi-squared  (top)  and  correlation  (bottom) 

values  (y-axis)  plotted  against  the  mean  satellite  elevation  angle  (x-axis) 

for  all  39  WV-2  images 
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Coastal  &  Green  Chi-Squared 


Coastal  &  Green  Correlation 


Mean  Satellite  Elevation  Angle 

Figure  37.  Coastal  &  Green:  Chi-squared  (top)  and  correlation  (bottom)  values 
(y-axis)  plotted  against  the  mean  satellite  elevation  angle  (x-axis)  for 

all  39  WV-2  images 
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Blue  &  Green  Chi-Squared 


Blue  &  Green  Correlation 


10%  1 - 1 - . - T - , - , - , - 7 - , - 7 - 7 - 7 - 7 - 

24.9  29.8  382  46.7  57.6  702  77.8  69.0  56.6  46.5  38.7  32.5  27.4 

Mean  Satellite  Elevation  Angle 


Figure  38.  Blue  &  Green:  Chi-squared  (top)  and  correlation  (bottom)  values 
(y-axis)  plotted  against  the  mean  satellite  elevation  angle  (x-axis)  for 

all  39  WV-2  images 
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4.40E+06 


Green  &  Yellow  Chi-Squared 
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Mean  Satellite  Elevation  Angle 


Green  &  Yellow  Correlation 


Figure  39.  Green  &  Yellow:  Chi-squared  (top)  and  correlation  (bottom)  values 
(y-axis)  plotted  against  the  mean  satellite  elevation  angle  (x-axis)  for 

all  39  WV-2  images 


Although  exact  chi-squared  and  correlation  values  depend  on  the  cutoff  values 
used  in  the  2D  scatter  plots  for  each  band  combination  (refer  to  Figure  23),  the  trend  was 
the  same  after  multiple  runs.  For  all  four  band  combinations,  the  images  acquired  closer 
to  nadir  had  the  lowest  chi-squared  values  and  the  highest  correlations.  The  optimal 
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values  were  found  when  running  the  Blue  &  Green  code,  with  a  correlation  value  around 
71%.  The  most  erratic  results  came  from  the  Green  &  Yellow  code.  This  may  be 
because  the  Yellow  band  did  not  penetrate  as  deeply  into  the  water. 

Low  chi-squared  values  are  also  seen  from  Image  IDs  4050  to  4100.  This  is  most 
likely  due  to  the  fact  that  clouds  were  entering  the  scene,  being  masked  and,  therefore, 
leaving  very  little  water  from  which  to  collect  spectral  infonnation  for  depth  derivation. 
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V.  CONCLUSIONS 


A.  SUMMARY  OF  OBSERVATIONS  AND  ANALYSIS 

Thirty-nine  WV-2  images  collected  at  multiple  angles  over  a  coastal,  shallow 
water  environment  were  used  to  analyze  the  effect  of  varying  view  angle  on  bathymetry 
derivation.  Following  initial  data  preparation  steps,  the  entire  dataset  was  analyzed  to 
determine  what  effect,  if  any,  acquisition  angle  has  on  nearshore  depth  retrieval.  The 
data  processing  and  analysis  consisted  of  steps  to  remove  sun  glint,  a  band  ratio  method 
to  determine  relative  depth,  and  a  regression  to  find  derived  depth. 

Accuracy  of  depth  retrieval  did,  in  fact,  demonstrate  an  association  with  image 
acquisition  angle.  Images  acquired  at  more  off-nadir  view  angles  proved  to  have  lower 
correlation  values  to  the  actual  “true”  bathymetry  data.  This  was  shown  by  the  increase 
in  correlation  values  for  more  nadir  images  within  the  dataset,  followed  by  a  decrease  as 
images  became  more  backward-looking.  Likewise,  chi-squared  values  decreased  as  the 
code  approached  the  nadir-looking  images,  and  increased  again  after  that. 

Four  band  combinations  were  compared  in  an  effort  to  detennine  which  WV-2 
bands  might  be  best  for  determining  depth  in  shallow,  coastal  environments.  These 
included:  Coastal  &  Blue,  Coastal  &  Green,  Blue  &  Green,  and  Green  &  Yellow.  Of  the 
four,  the  Blue  &  Green  band  combination  perfonned  best,  with  the  highest  correlation 
value  of  approximately  71%  between  “true”  and  derived  depth. 

The  results  of  this  research  showed  that  more  accurate  bathymetric  depths  will  be 
derived  from  images  that  have  been  acquired  closer  to  nadir  viewing  geometry. 
Additionally,  when  using  imagery  from  the  WV-2  sensor,  a  combination  of  the  Blue  and 
Green  bands  will  penetrate  into  the  water  in  such  a  way  that  a  higher  accuracy  of  depth 
information  will  be  obtained. 

B.  RECOMMENDATIONS  AND  FUTURE  WORK 

There  are  a  number  of  factors  that  may  have  negatively  affected  the  data  and 
analysis.  While  the  scene  over  Kailua  Bay  was  chosen  because  it  was  the  least  cloudy 
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beach  out  of  the  39  images  in  the  dataset,  it  still  had  cloud  cover  (especially  between 
Image  IDs  4050  and  4100).  This  led  to  fewer  data  points  that  could  be  analyzed  for  those 
scenes.  It  would  also  have  been  preferable  to  have  used  a  more  recent  dataset  for  “true” 
bathymetry.  This  particular  collection  is  thought  to  have  been  acquired  between  2002 
and  2005,  and  in  an  environment  as  dynamic  as  a  coastal  shoreline,  bathymetry  could 
have  changed  significantly.  In  the  future,  it  would  be  advantageous  to  find  a  recent 
bathymetric  dataset,  and  then  acquire  WV-2  data  over  that  geographic  location. 

There  are  a  number  of  possibilities  for  improvement  and  growth  in  future  research 
related  to  this  particular  topic.  It  would  be  beneficial  to  apply  the  analysis  approach  to 
the  entire  dataset  scene,  rather  than  just  the  image  chips.  It  would  also  be  interesting  to 
use  a  dataset  with  a  wider  range  of  satellite  elevation  angles  (more  off-nadir  to  more 
nadir  and  back).  It  might  also  be  fascinating  to  apply  this  method  to  nearshore  areas  with 
black  sand  to  see  whether  or  not  sand  color  affects  the  results  of  the  algorithm. 


58 


LIST  OF  REFERENCES 


Avery,  T.  E.,  G.  L.  Berlin  (1992),  Fundamentals  of  Remote  Sensing  and  Airphoto 

Interpretation,  Macmillan  Publishing  Company,  Maxwell  Macmillan  Canada, 
Maxwell  Macmillan  International,  New  York,  Toronto,  Oxford,  Singapore, 
Sydney. 

Camacho,  M.  A.  (2006),  Depth  Analysis  of  Midway  Atoll  Using  QuickBird  Multi- 
Spectral  Imaging  over  Variable  Substrates,  M.S.  thesis,  Naval  Postgraduate 
School,  Monterey,  CA. 

Clark,  R.  E.  (2005),  Naval  Satellite  Bathymetry:  A  Performance  Assessment,  M.S.  thesis, 
Naval  Postgraduate  School,  Monterey,  CA. 

Dickey,  T.  D.,  G.  W.  Kattawar,  and  K.  J.  Voss  (2011),  Shedding  new  light  on  light  in  the 
ocean,  Physics  Today  (April),  44. 

Eismann,  M.  T.  (2012),  Hyperspectral  Remote  Sensing,  SPIE  Press,  Bellingham,  WA. 

Goodman,  J.  A.,  Z.  Lee,  and  S.  L.  Ustin  (2008),  Influence  of  Atmospheric  and  Sea- 

Surface  Corrections  on  Retrieval  of  Bottom  Depth  and  Reflectance  Using  a  Semi- 
Analytical  Model:  A  Case  Study  in  Kaneohe  Bay,  Hawaii,  Applied  Optics,  47,  FI. 

Hedley,  J.  D.,  A.  R.  Harbome,  and  P.  J.  Mumby  (2005),  Simple  and  robust  removal  of 

sun  glint  for  mapping  shallow- water  benthos,  Int.  J.  Remote  Sens.,  26(10),  2107. 

Hochberg,  E.  J.,  S.  Andrefouet,  and  M.  R.  Tyler  (2003),  Sea  Surface  Correction  of  High 
Spatial  Resolution  Ikonos  Images  to  Improve  Bottom  Mapping  in  Near-Shore 
Environments,  IEEE  Trans.  Geosci.  Remote  Sens.,  41,  1724. 

Isoun,  E.,  C.  H.  Fletcher,  N.  Frazer,  and  J.  Gradie  (2003),  Multi-spectral  mapping  of  reef 
bathymetry  and  coral  cover;  Kailua  Bay,  Hawaii,  Coral  Reefs,  22,  68, 
doi:  10. 1007/s00338-003-0287-4. 

Jensen,  J.  R.  (2007),  Remote  Sensing  of  the  Environment:  An  Earth  Resource 
Perspective,  vol.  2,  Prentice  Hah,  Upper  Saddle  River,  NJ. 

Kay,  S.,  J.  D.  Hedley,  and  S.  Lavender  (2009),  Sun  Glint  Correction  of  High  and  Low 
Spatial  Resolution  Images  of  Aquatic  Scenes:  a  Review  of  Methods  for  Visible 
and  Near-Infrared  Wavelengths,  Remote  Sensing,  7(4),  697. 

Longbotham,  N.,  C.  Bleiler,  C.  Chaapel,  C.  Padwick,  W.  Emery,  and  F.  Pacifici  (2011), 
Spatial  Classification  of  WorldView-2  Multi-angle  Sequence,  paper  presented  at 
Joint  Urban  Remote  Sensing  Event,  Munich,  Germany. 


59 


Longbotham,  N.,  C.  Bleiler,  C.  Chaapel,  C.  Padwick,  W.  Emery,  and  F.  Pacifici  (2011), 
Spectral  Classification  of  WorldView-2  Multi-angle  Sequence,  paper  presented  at 
Joint  Urban  Remote  Sensing  Event,  Munich,  Germany. 

Longbotham,  N.,  C.  Chaapel,  L.  Bleiler,  C.  Padwick,  W.  J.  Emery,  and  F.  Pacifici  (2011), 
Very  High  Resolution  Multiangle  Urban  Classification  Analysis,  IEEE 
TRANSACTIONS  ON  GEOSCIENCE  AND  REMOTE  SENSING(99),  1 , 
doi:10.1109/TGRS.201 1.2165548. 

Loomis,  M.  J.  (2009),  Depth  Derivation  from  the  WorldView-2  Satellite  Using 

Hyperspectral  Imagery,  M.S.  thesis,  Naval  Postgraduate  School,  Monterey,  CA. 

Lyzenga,  D.  R.  (1978),  Passive  remote  sensing  techniques  for  mapping  water  depth  and 
bottom  features,  Applied  Optics,  77(3),  379,  doi:10.1364/AO.  17.000379. 

Lyzenga,  D.  R.,  N.  P.  Malinas,  and  F.  J.  Tanis  (2006),  Multispectral  Bathymetry  Using  a 
Simple  Physically  Based  Algorithm,  IEEE  Trans.  Geosci.  Remote  Sens.,  44, 
2251. 

Madden,  C.  K.  (2011),  Contributions  to  Remote  Sensing  of  Shallow  Water  Depth  with 
the  Worldview-2  Yellow  Band,  M.S.  thesis,  Naval  Postgraduate  School, 
Monterey,  CA. 

Mancini,  S.  (2012),  Automating  Nearshore  Bathymetry  Extraction  from  Wave  Motion  in 
Satellite  Optical  Imagery,  M.S.  thesis,  Naval  Postgraduate  School,  Monterey,  CA. 

Marchisio,  G.,  F.  Pacifici,  and  C.  Padwick  (2011),  Coastal  applications  of  WorldView-2 
high  resolution  multi-spectral  imagery,  paper  presented  at  American  Society  for 
Photogrammetry  and  Remote  Sensing,  Milwaukee,  WI  USA. 

Maritorena,  S.,  A.  Morel,  and  B.  Gentili  (1994),  Diffuse  reflectance  of  oceanic  shallow 
waters:  Influence  of  water  depth  and  bottom  albedo,  Limnol.  Oceanogr.,  39(1), 
1689. 

Martin,  S.  (2004),  An  Introduction  to  Ocean  Remote  Sensing,  Cambridge  University 
Press,  Cambridge,  UK. 

McCarthy,  B.  L.  (2010),  Coastal  Bathymetry  Using  8-Color  Multispectral  Satellite 

Observation  of  Wave  Motion,  M.S.  thesis,  Naval  Postgraduate  School,  Monterey, 
CA. 

McConnon,  C.  L.  (2010),  High  Spatial  Resolution  Bidirectional  Reflectance  Retrieval 
Using  Satellite  Data,  M.S.  thesis,  Naval  Postgraduate  School,  Monterey,  CA. 

Mobley,  C.  D.  (1994),  Light  and  Water:  Radiative  Transfer  in  Natural  Waters,  Academic 
Press,  San  Diego,  CA  USA. 


60 


Myrick,  K.  B.  (2011),  Coastal  Bathymetry  Using  Satellite  Observation  in  Support  of 
Intelligence  Preparation  of  the  Environment,  M.S.  thesis,  Naval  Postgraduate 
School,  Monterey,  CA. 

Olsen,  R.  C.  (2007),  Remote  Sensing  from  Air  and  Space,  SPIE  Press,  Bellingham, 
Washington. 

Pacifici,  F.,  J.  Chanussot,  and  Q.  Du,  2011:  2011  GRSS  Data  Fusion  Contest:  Exploiting 
WorldView -2  Multi- Angular  Acquisitions. 

Pacifici,  F.,  K.  Navular,  2011:  Very  High  Spatial  Resolution  Imagery  for  Urban 

Application.  [Available  online  at  http://www.earthzine.org/201 1/06/27/very-high- 
spatial-resolution-imagery-for-urban-applications/]. 

Research  Systems,  I.  (2004),  ENVI  Tutorials,  Boulder,  CO:  Research  Systems,  Inc. 

Sanford,  T.  B.,  K.  A.  Kelly,  and  D.  M.  Farmer  (2011),  Sensing  the  ocean,  Physics  Today 
(February),  24. 

Smith,  R.  C.,  K.  S.  Baker  (1981),  Optical  properties  of  the  clearest  natural  waters, 
Applied  Optics,  20(2),  111 . 

Stumpf,  R.  P.,  K.  Holderied,  and  M.  Sinclair  (2003),  Detennination  of  water  depth  with 
high-resolution  satellite  imagery  over  variable  bottom  types,  Limnol.  Oceanogr., 
48(  1,  part  2),  547. 

Sturm,  B.  (1981),  The  atmospheric  correction  of  remotely  sensed  data  and  the 

quantitative  determination  of  suspended  matter  in  marine  water  surface  layers  in 
Remote  Sensing  in  Meteorology,  Oceanography  and  Hydrology,  edited  by  A.  P. 
Cracknell,  ed.  (Ellis  Horwood,  Chichester,  UK,  1980),  pp.  163-197. 

Thomas,  G.  E.,  K.  Stamnes  (1999),  Radiative  Transfer  in  the  Atmosphere  and  Ocean, 
Cambridge  University  Press,  Cambridge,  UK. 

Updike,  T.,  C.  Comp  (2010),  Radiometric  Use  of  WorldView-2  Imagery,  technical  note, 
rev.  1.0,  DigitalGlobe,  Inc.,  Fongmont,  CO. 


61 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


62 


APPENDIX:  CHARTS  AND  GRAPHS 


A.  EXAMPLE  WORLDVIEW-2  METADATA  FILE  (FOR  IMAGE  1010) 


version  =  "23.6"; 

generationTime  =  2011-08-01T01 : 59 : 55 . 000000Z; 

productOrderld  =  "052564471010  01  P001"; 

productCatalogld  =  "202001008362B200" ; 

imageDescriptor  =  "BasiclB"; 

bandld  =  "Multi"; 

panSharpenAlgorithm  =  "None"; 

numRows  =  3072; 

numColumns  =  9216; 

productLevel  =  "LV1B"; 

productType  =  "Basic"; 

numberOf Looks  =  1; 

radiometricLevel  =  "Corrected"; 

bitsPerPixel  =  16; 

compressionType  =  "JPEG2000"; 

jpegProf ileName  =  "nga  npje_pan  nl"; 

BEGINJGROUP  =  BAND_C 

-157.86756401; 

21.36508648; 

275.36; 

-157 .49196245; 

21.69660890; 

9.00; 

-157.49051787; 

21.56876101; 

11 .00; 

-157 . 87041130; 

21 . 22582540; 

15 .00; 

absCalFactor  =  9 . 2 95654e-03 ; 

ef f ectiveBandwidth  =  4 . 730000e-02 ; 

TDILevel  =  24; 

END_GROUP  =  BAND_C 
BEGIN  GROUP  =  BAND  B 


ULLon 

ULLat 

ULHAE 

URLon 

URLat 

URHAE 

LRLon 

LRLat 

LRHAE 

LLLon 

LLLat 

LLHAE 


ULLon 

=  -157. 

86756401; 

ULLat 

=  21 . 

36508648; 

ULHAE 

=  275 

.36; 

URLon 

=  -157. 

49196245; 

URLat 

=  21 . 

69660890; 

URHAE 

9 

o 

o 

LRLon 

=  -157. 

49051787; 

LRLat 

=  21 . 

56876101; 

LRHAE 

=  11 

o 

o 

LLLon 

=  -157. 

87041130; 

LLLat 

=  21 . 

22582540; 
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LLHAE  =  15.00; 

absCalFactor  =  1 . 783568e-02 ; 

ef f ectiveBandwidth  =  5 . 430000e-02 ; 

TDILevel  =  10; 

END_GROUP  =  BAND_B 
BEGIN_GROUP  =  BAND_G 

ULLon  =  -157.86756401; 

ULLat  =  21.36508648; 

ULHAE  =  275.36; 

URLon  =  -157.49196245; 

URLat  =  21 . 69660890; 

URHAE  =  9.00; 

LRLon  =  -157.49051787; 

LRLat  =  21.56876101; 

LRHAE  =  11.00; 

LLLon  =  -157.87041130; 

LLLat  =  21.22582540; 

LLHAE  =  15.00; 

absCalFactor  =  1 . 3641 97e-02 ; 
eff ectiveBandwidth  =  6 . 300000e-02 ; 
TDILevel  =  10; 

END_GROUP  =  BAND_G 
BEGIN_GROUP  =  BAND_Y 

ULLon  =  -157.86756401; 

ULLat  =  21.36508648; 

ULHAE  =  275.36; 

URLon  =  -157.49196245; 

URLat  =  21 . 69660890; 

URHAE  =  9.00; 

LRLon  =  -157.49051787; 

LRLat  =  21.56876101; 

LRHAE  =  11.00; 

LLLon  =  -157.87041130; 

LLLat  =  21.22582540; 

LLHAE  =  15.00; 

absCalFactor  =  6 . 81 07 18e-03 ; 
eff ectiveBandwidth  =  3 . 740000e-02 ; 
TDILevel  =  18; 

END_GROUP  =  BAND_Y 
BEGIN_GROUP  =  BAND_R 

ULLon  =  -157.86756401; 

ULLat  =  21.36508648; 

ULHAE  =  275.36; 

URLon  =  -157.49196245; 

URLat  =  21 . 69660890; 

URHAE  =  9.00; 

LRLon  =  -157.49051787; 

LRLat  =  21.56876101; 

LRHAE  =  11.00; 

LLLon  =  -157.87041130; 
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LLLat  =  21.22582540; 

LLHAE  =  15.00; 

absCalFactor  =  1 . 851735e-02 ; 

ef f ectiveBandwidth  =  5 . 740000e-02 ; 

TDILevel  =  6; 

END_GROUP  =  BAND_R 
BEGIN_GROUP  =  BAND_RE 

ULLon  =  -157.86756401; 

ULLat  =  21.36508648; 

ULHAE  =  275.36; 

URLon  =  -157.49196245; 

URLat  =  21 . 69660890; 

URHAE  =  9.00; 

LRLon  =  -157.49051787; 

LRLat  =  21.56876101; 

LRHAE  =  11.00; 

LLLon  =  -157.87041130; 

LLLat  =  21.22582540; 

LLHAE  =  15.00; 

absCalFactor  =  6 . 0631 45e-03 ; 
eff ectiveBandwidth  =  3 . 930000e-02 ; 
TDILevel  =  18; 

END_GROUP  =  BAND_RE 
BEGIN_GROUP  =  BAND_N 

ULLon  =  -157.86756401; 

ULLat  =  21.36508648; 

ULHAE  =  275.36; 

URLon  =  -157.49196245; 

URLat  =  21 . 69660890; 

URHAE  =  9.00; 

LRLon  =  -157.49051787; 

LRLat  =  21.56876101; 

LRHAE  =  11.00; 

LLLon  =  -157.87041130; 

LLLat  =  21.22582540; 

LLHAE  =  15.00; 

absCalFactor  =  2 . 050828e-02 ; 
eff ectiveBandwidth  =  9 . 890000e-02 ; 
TDILevel  =  6; 

END_GROUP  =  BAND_N 
BEGIN_GROUP  =  BAND_N2 

ULLon  =  -157.86756401; 

ULLat  =  21.36508648; 

ULHAE  =  275.36; 

URLon  =  -157.49196245; 

URLat  =  21 . 69660890; 

URHAE  =  9.00; 

LRLon  =  -157.49051787; 

LRLat  =  21.56876101; 

LRHAE  =  11.00; 
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LLLon  =  -157.87041130; 

LLLat  =  21.22582540; 

LLHAE  =  15.00; 

absCalFactor  =  9 . 042234e-03 ; 
ef f ectiveBandwidth  =  9 . 960000e-02 ; 
TDILevel  =  24; 


END_GROUP  =  BAND_N2 
outputFormat  =  "NITF2 1NCDRD" ; 

BEGIN_GROUP  =  IMAGE_1 
satld  =  "WV02"; 
mode  =  "FullSwath"; 
scanDirection  =  "Forward"; 

Catld  =  " 103001 00  0DC522  00 " ; 

TLCTime  =  2011-07-30T21 : 22 : 48 . 765050Z; 
numTLC  =  2 ; 

TLCList  =  ( 

(0,  0.000000), 

(2055,  0.411000)  ) ; 

=  2011-07-30T21 : 22 : 48 . 765050Z; 


firstLineTime 


avgLineRate  =  5000.00; 
exposureDuration  =  0.0002; 
minCollectedRowGSD  =  7.002; 

maxCollectedRowGSD  =  7.136; 

meanCollectedRowGSD  =  7.069; 

minCollectedColGSD  =  6.088; 

maxCollectedColGSD  =  6.168; 


meanCollectedColGSD  = 

6.128; 

meanCollectedGSD  = 

6.582; 

meanProductRowGSD 

= 

7 .110; 

meanProductColGSD 

= 

6.162; 

meanProductGSD  = 

6. 

619; 

rowUncertainty  = 

294 

o 

LO 

colUncertainty  = 

145 

.97; 

minSunAz  = 
maxSunAz  = 
meanSunAz  = 
minSunEl  = 
maxSunEl  = 
meanSunEl  = 
minSatAz  = 
maxSatAz 


; 

4; 


meanSatAz 


minSatEl  = 
maxSatEl  = 


meanSatEl 


96. 

96. 

96.2; 

72.3; 

72.3; 

72.3; 
18. 

18. 

1J 
24  . 

25. 

24 . 9; 


9; 

0; 


7; 

9; 


minlnTrackViewAngle  =  52 

maxInTrackViewAngle  =  52 

meanlnTrackViewAngle  =  52.8; 

minCrossTrackViewAngle  =  14, 

maxCrossTrackViewAngle  =  14, 


3; 

4; 
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meanCrossTrackViewAngle  =  14.4; 

minOf fNadirViewAngle  =  54.1; 

=  54.1; 

=  54.1; 


cloudCover  =  0.311; 
resamplingKernel  =  "CC"; 
positionKnowledgeSrc  =  "R"; 
attitudeKnowledgeSrc  =  "R"; 
revNumber  =  9481; 

END_GROUP  =  IMAGE_1 
END; 


maxOf fNadirViewAngle 
meanOf fNadirViewAngle 
PNIIRS  =  1.3; 
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B.  INFORMATION  ABOUT  THE  39  WORLDVIEW-2  IMAGES  -  IMAGE  ID,  ACQUISITION  TIME,  SENSOR 
ELEVATION,  SENSOR  AZIMUTH,  OFF  NADIR  VIEW  ANGLE,  AZIMUTH,  AND  ZENITH 


Image 

ID 

firstLineTime 

(Z) 

meanSatEl 

meanSatAz 

meanOffN  adirViewAngle 

Azimuth 

Zenith 

1010 

21:22:49 

24.9 

18.8 

54.1 

-161.2 

125.9 

1020 

21:22:58 

26.5 

19.2 

53.1 

-160.8 

126.9 

1030 

21:23:07 

28.1 

19.6 

52.0 

-160.4 

128.0 

1040 

21:23:16 

29.8 

20.1 

50.8 

-159.9 

129.2 

1050 

21:23:25 

31.7 

20.6 

49.5 

-159.4 

130.5 

1060 

21:23:34 

33.7 

21.2 

48.0 

-158.8 

132.0 

1080 

21:23:53 

38.2 

22.6 

44.6 

-157.4 

135.4 

1090 

21:24:02 

40.8 

23.5 

42.6 

-156.5 

137.4 

1100 

21:24:12 

43.7 

24.6 

40.3 

-155.4 

139.7 

2010 

21:24:21 

46.7 

26.0 

37.8 

-154.0 

142.2 

2020 

21:24:31 

50.1 

27.6 

35.1 

-152.4 

144.9 

2030 

21:24:41 

53.6 

29.7 

32.1 

-150.3 

147.9 

2040 

21:24:51 

57.6 

32.4 

28.7 

-147.6 

151.3 

2050 

21:25:01 

61.6 

36.1 

25.2 

-143.9 

154.8 

2060 

21:25:12 

66.0 

41.5 

21.4 

-138.5 

158.6 

2070 

21:25:22 

70.2 

49.2 

17.4 

-130.8 

162.3 

2080 

21:25:32 

74.0 

61.2 

14.3 

-118.8 

165.7 

2090 

21:25:42 

76.9 

80.1 

11.7 

-99.9 

168.3 

2100 

21:25:52 

77.8 

105.5 

10.9 

-74.5 

169.1 

3010 

21:26:02 

76.3 

129.8 

12.1 

-50.2 

167.9 

3020 

21:26:12 

73.1 

146.9 

15.0 

-33.1 

165.0 

3030 

21:26:23 

69.0 

157.9 

18.6 

-22.1 

161.4 

68 


firstLineTime 


meanSatEl 

meanSatAz 

21:26:33 

64.8 

164.7 

21:26:43 

60.6 

169.4 

21:26:53 

56.6 

172.8 

21:27:02 

52.9 

175.3 

21:27:12 

49.6 

177.2 

21:27:21 

46.5 

178.7 

21:27:30 

43.7 

179.9 

21:27:39 

41.1 

180.9 

21:27:47 

38.7 

181.7 

21:27:56 

36.5 

182.5 

21:28:04 

34.4 

183.1 

21:28:13 

32.5 

183.7 

21:28:21 

30.7 

184.2 

21:28:29 

29.0 

184.6 

21:28:38 

27.4 

185.0 

21:28:46 

25.9 

185.4 

21:28:54 

24.5 

185.7 

c. 


IDL  CODE  FOR  MOST  HIGHLY  CORRELATED  BAND  COMBINATION  (BLUE  &  GREEN) 


;Final  code  for  Krista  Lee's  September  2012  thesis 
;Code  written  by:  Krista  Lee,  Angie  Kim,  Professor  Olsen 

r 

;Blue  vs.  Green 

r 

;The  purpose  of  this  code  is  to  take  39  image  chips  over  Kailua  Bay,  Oahu,  HI  from  30  July  2011  and  apply  a 
deglinting  process, 

; along  with  a  band  ratio  method  to  determine  water  depth.  The  images,  acquired  at  different  angles,  have 
alreaad  been  calibrated 

;for  radiance,  atmospherically  corrected,  and  registered  to  the  most  nadir  of  the  images.  This  code  will 
then  apply  a  simple 

;deglinting  technique  to  remove  the  effects  of  sun  glint,  as  well  as  use  Stumpf  et  al . ' s  band  ratio  method 
for  bathymetry  derivation. 

; These  depth  values  will  be  compared  to  actual  lidar  bathymetry  data  from  the  University  of  Hawaii. 

•  •k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

r 

; LOAD  COLOR  TABLE 

loadct,  13  ;  Loads  color  table  --  Rainbow 

tvlct,  r,  g,  b,  /get 
r  (255)  =  255 
g  (255)  =  255 
b  (255)  =  255 
tvlct,  r,  g,  b 

device,  decomposed  =  0,  retain  =  1 

•  -k'k-k'k-k'k'k'k-k'k'k'k-k'k'k'k'k'k-k'k-k'k-k'k'k'k'k-k-k'k-k'k-k'k'k'k'k'k-k'k-k'k'k'k-k'k'k'k-k'k-k'k'k'k'k'k 

r 

; DEFINE  DIRECTORY 

dir  =  'O:\Krista  Lee\Thesis\IDL\FINAL ' 

cd,  dir 

•  'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

r 

; DEFINE  FILES 

files  =  FILE_SEARCH ( ' *  mosaic  wgs84  rad  resize*chip . dat ' ) 

nfiles  =  n_elements ( files ) 

print,  string (nfiles )  +  '  files  found:  ' 
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print,  files 


fmask  =  '2100  mask  b.dat'  ;Created  specifically  for  this  image  chip  in  ENVI 

;Size  of  images 

samples  =  995L 
lines  =  999L 
bands  =  8 

; FOR  LOOP  TO  REPEAT  PROCESS  FOR  EACH  IMAGE 

for  file  number  =  0,  nfiles  -  1  do  begin  ; Loops  through  each  elements  in  file  array 

file  =  files (file  number) 

out_dir  =  'Image_'  +  strmid (file,  0,  4) 

spawn,  'mkdir  '  +  out  dir  ;This  creates  the  folder  for  the  WV-2  image 
cd,  out_dir 

spawn,  'mkdir  b  g'  /This  creates  the  next  folder  down  for  the  band  combinations 

cd,  'b_g',  CURRENT  =  old^dir 

help,  old_dir 

outdir  =  old_dir  +  ' \b_g\ ' 

cd,  dir  ; Go  back  to  the  main  directory  where  all  the  datafiles  are,  etc. 

/DEFINE  OUTPUT  FILES 

/Use  input  filename  to  create  output  filename 

filel9  =  files (file  number)  /String  containing  name  of  the  current  file  (filel9  =  current  file) 

shortf ilename  =  strmid (filel9,  0,  4) 

band  ratio  vs  true  depth  =  shortf ilename  +  '  band  ratio  vs  true  depth.tif 
regression  file  =  shortf ilename  +  '  regression.tif 
final  depth  file  =  shortf ilename  +  '  final  depth.tif' 

der  vs  true  scatter  plot  file  =  shortf ilename  +  '  der  vs  true  scatter_plot.tif 

•  -k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k-k'k-k'k-k'k'k'k 

r 

/OPEN  LAND  MASK 

mask  =  by tarr (samples,  lines) 
openr,  1,  fmask 
forrd,  1,  mask 
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close,  1 


land  mask  =  mask 

mask  =  reform (mask,  samples,  lines) 
mask  =  rebin (mask,  samples,  lines,  bands) 

datal  =  fltarr (samples,  lines,  bands) 
openr,  1,  filel9 
forrd,  1,  datal 

close,  1 

; DISPLAY  LAND-MASKED  IMAGE 

window,  1,  xsize  =  samples,  ysize  =  lines,  title  =  '  (1)  Water  Only  (land  masked)  -  '+filel9 
datal  =  datal  *  mask 

tvscl,  datal (*,*,  [4,2, 1] )  ,  true  =  3,  order  =  1 

; MASK  ADDITIONAL  PROBLEM  PIXELS  (GLINT,  CLOUDS,  ETC.) 

;User  defines  water  region 

window,  2,  title  =  ' (2)  Scatter  Plot  -  Blue  vs.  NIR-1  -  '  +  filel9 

ti  =  datal (*,*,1)  ;Blue 
tj  =  datal (*, *, 6)  ;NIR1 

order  =  sort(ti)  &  ti  =  ti (order)  &  tj  =  tj (order)  ;Not  strictly  necessary,  makes  oplot  easier 
nele  =  n_elements (order)  &  index  =  10*  lindgen (nele/10)  ;  Plots  every  10th  point 

plot,  ti (  index),  tj (index) ,  psym  =  3,  title  =  file,  xtitle  =  'Blue',  ytitle  =  'NIR-1' 

print,  'Click  the  upper-right-hand  corner  of  the  region  that  defines  water  --  for  masking  purposes' 

cursor,  x,  y,  /down,  /data 
delta  =  0.1 

xx  =  [0,  x,  x,  0,  0]  ;For  lines  that  create  box  when  clicked 

yy  =  [  0,  0,  y  ,  y,  0] 

plots,  xx,  yy 

glint_clouds_etc  =  where  (datal (*,*,  1 )  gt  x  or  datal  (*,  *,  6)  gt  y) 
gee  =  glint_clouds_etc 

ngee  =  where (datal (*,*, 1 )  It  x  and  datal (*,*,6)  It  y)  ; ngee :  NOT  glint,  clouds,  etc. 
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help,  gee 

datal=  reform (datal,  samples*lines,  bands) 
datal (gee,  * )  =0 

datal  =  reform (datal ,  samples,  lines,  bands) 

gce_mask  =  by tarr  ( samples ,  lines)  &  gce_mask(*)  =  1  &  gce_mask (gee)  =  0 
land  gee  mask  =  land  mask*gce  mask 

window,  3,  xsize  =  samples,  ysize  =  lines,  title  =  ' (3)  Water  Only  (land  and  clouds  masked)  -  '  +  filel9 
tvscl,  datal (*,*,  [4,2, 1] ) ,  true  =  3,  order  =  1 

*-k'k'k'k-k'k-k'k-k'k'k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k-k'k-k'k-k'k'k'k-k'k'k'k'k'k-k'k-k'k-k'k-k'k-k'k 

r 

; DEGLINT  USING  HEDLEY  APPROACH  (UTILIZES  SLOPE) 

; Shows  scatter  plot,  best  fit  line,  slop,  and  line  along  the  length  of  the  range 

windo’  ,  4,  title  =  ' (4)  Scatter  plot  and  fit  for  deglint  (NIR-1  vs.  Blue)  -  '  +  filel9 
ti  =  datal  (*,*,1)  ,-Blue 

tj  =  datal (*, *, 6)  ;NIR1 

; Check  for  zeroes 

zeroes  =  where (ti  le  0) 
nz  =  where (ti  gt  0) 

ti  =  ti (nz) 
tj  =  tj  (nz) 

order  =  sort(ti)  &  ti  =ti (order)  &  tj  =  tj (order)  ;Not  strictly  necessary,  makes  oplot  easier 

nele  =  n_elements (  order)  &  index  =  10  *  lindgen (  nele/10)  ; Plots  every  10th  point 

plol  ,  tj  (index) ,  ti(index),  psym  =  3,  /ynozero,  xrange  =  [0,  2.0],  xtitle  =  'NIR-1',  ytitle  =  'Blue' 
degree=  1 

coeff_b  =  poly_£  (tj,  ti,  degree,  chisq=chisql ,  yfit  =  yfitb) 

oplot,  t j ,  yfitb 

slope_b  =  coeff_b(l) 

xyouts,  0.2,  .2,  string ( slope_b) ,  /normal 

water_irl_min  =  min(tj) 
water  irl  max  =  max(tj) 
xmin  =  min(ti) 

plots,  [water  irl  min,  water  irl  max] ,  [xmin,  xmin] 

datal (*,*,1)  =  datal (*,*,1)  -  slope_b  *  (datal ( * , * , 6)  -  water__irl_min)  ;Blue  and  NIR-1  bands 

window,  5,  title  =  ' (5)  Scatter  plot  and  fit  for  deglint  (NIR-1  vs.  Green)  - 
ti  =  datal (*,*,2)  ; Green 


76 


+  f ilel9 


;  NIR1 


t j  =  datal ( * ,  * ,  6) 

; Check  for  zeroes 

zeroes  =  where (ti  le  0) 
nz  =  where (ti  gt  0) 

ti  =  ti (nz) 
tj  =  tj  (nz) 

order  =  sort(ti)  &  ti  =  ti (order)  &  tj  =  tj (order)  ;Not  strictly  necessary,  makes  oplot  easier 

nele  =  n_elements (  order)  &  index  =  10  *  lindgen (nele/10)  ; Plots  every  10th  point 

plot,  tj (index) ,  ti(index),  psym  =  3,  /ynozero,  xrange  =  [0,  2.0],  xtitle  =  'NIR-1',  ytitle  =  'Green' 
degree  =  1 

coeff_g  =  poly_f  (  t j ,  ti,  degree,  chisq=chisql ,  yfit  =  yfitg) 

oplot,  t j ,  yfitg 

slope_g  =  coeff_g(l) 

xyouts,  0.2,  .2,  string ( slope_g) ,  /normal 

water_irl_min  =  min(tj) 
water  irl  max  =  max(tj) 
xmin  =  min(ti) 

plots,  [water  irl  min,  water  irl  max] ,  [xmin,  xmin] 

datal (*,*,2)  =  datal (*,*,2)  -  slope_g  *  (datal ( * ,  * ,  6)  -  water_irl_min)  ;Green  and  NIR-1  bands 

********************************************************* 

r 

; DEGLINT  IMAGERY  BY  SUBTRACTING  IR  WAVELENGTH  REFLECTANCE 
imlb  =  datal 

•  'k'k'k'k-k'k-k'k'k'k'k'k-k'k-k'k-k'k'k'k'k'k'k'k'k'k'k'k-k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

r 

; OBTAIN  RELATIVE  BATHYMETRY  BY  RATIO  OF  BANDS  METHOD 

;Use  expression:  In (1000*bl) /In (1000*b2) ,  where  bl  =  Blue  band,  b2  =  Green  band 
indexl  =  where (imlb  It  1)  &  imlb(indexl)  =  1 

nfac  =  100.  ;Camacho  used  1000. 

reldepth  =  alog (nfac* imlb (*,*, 1 )) /alog (nfac* imlb (*,*, 2 ) ) 

window,  6,  xsize  =  samples,  ysize  =  lines,  title  =  ' (6)  Blue/Green  Band  Ratio  -  Bathymetry  -  '  +  filel9 
rl  =  -1.1  ;min (truedepth) 

r2  =  -1.0  ;max (truedepth) 

tv,  bytscl (reldepth*land  gee  mask* (-1) ,  min=rl,  max=r2),  order  =  1 

bar  =  bindgen(256) 

bar  =  reform(bar,  1,  256) 

bar  =  rebin (bar,  40,  256) 
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xO  =  65  ;xO,  yO :  define  location  of  color  bar  within  the  image 
yO  =  30 

xl  =  xO  +  40;  +  1  ; xl ,  yl :  define  location  of  white  box  drawn  around  color  bar 

yl  =  yO  +  255;  +  1 
tv,  bar,  xO,  yO 

; Displays  color  bar  on  the  image 

plo'  ,  bar,  /nodata,  /noerase,  color  =  0,  xrange  =  [0,  1],  yrange  =  [rl,  r2]  *  (-1),  $ 
xticks  =  1,  yticks  =  4,  yminor  =  1,  charsize  =  1.0,  xstyle  =  1,  ystyle  =  1,  /device,  $ 
pos  =  [xO-l,  y0-l,  xl,  yl],  xtickname  =  ['  '],  yticklen  =  -0.05 

xyouts,  xl+25,  yO+128,  'Blue/Green  Band  Ratio',  orient  =  90,  size  =  1.0,  color  =  0,  /device,  align  = 

0.5 


; Close  unnecessary  windows 

wdelete,  1  ;Closes  Water  Only  (land  masked)  window 
;  wdelete,  4  ; Closes  NIR-1  vs.  Green  scatter  plot 

;  wdelete,  5  ;Closes  NIR-2  vs.  Yellow  scatter  plot 

•  -k'k-k'k-k'k-k'k-k-k'k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k-k'k'k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k 

r 

; OPEN  GROUNDTRUTHED  WATER  DEPTH  IMAGE 

truedepth  file  =  ' KailuaRaster  DepthMeters  resize  WGS84  2mpixels  995x999.dat' 

truedepth  =  dblarr ( samples ,  lines) 
openr,  1,  truedepth_f ile 
forrd,  1,  truedepth 

close,  1 

truedepth  =  truedepth  *  (-1.0)  ;Make  depth  values  positive  to  match  reldepth 

rl  =  -18  ;min (truedepth) 
r2  =  0  ;max (truedepth) 

bad  =  where (truedepth  gt  rl  *  (-1))  ;Gets  rid  of  values  deeper  than  -20  m  (rl)  --  assume  these  are  bad 

data 

nodata_mask  =  by tarr ( samples ,  lines)  &  nodata_mask ( * )  =  1  &  nodata_mask (bad)  =  0 
land  gee  mask  =  land  gee  mask  *  nodata  mask 

; Place  all  masks  on  truedepth 
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(-1) 


truedepth  masked  =  truedepth  *  land  gee  mask  * 
truedepth  =  truedepth  masked 

goto,  skipgt  /Skip  ground  truth  image  --  comment  out  this  line  if  you  want  a  window 

window,  7,  xsize  =  samples,  ysize  =  lines,  title  =  ' (7)  True  Bathymetry  -  '  +  truedepth  file 

tv,  bytscl (truedepth,  min  =  rl,  max  =  r2),  order  =  1 

bar  =  bindgen(256)  &  bar  =  reform(bar,  1,  256)  &  bar  =  rebin (bar,  40,  256) 

xO  =  65  ;x0,  yO :  define  location  of  color  bar  within  the  image 

yO  =  30 

xl  =  xO  +  40;  +  1  ;xl,  yl :  define  location  of  white  box  drawn  around  color  bar 

yl  =  yO  +  255;  +  1 
tv,  bar,  xO,  yO 

; Displays  color  bar  on  the  image 

ploi  ,  bar,  /nodata,  /noerase,  color  =  0,  xrange  =  [0,  1],  yrange  =  [rl,  r2]*(-l),  $ 
xticks  =  1,  yticks  =  3,  yminor  =  1,  charsize  =  1.0,  xstyle  =  1,  ystyle  =  1,  /device,  $ 
pos  =  [xO-l,  yO-l,  xl,  yl],  xtickname  =  ['  ,  yticklen  =  -0.05 

xyouts,  xl+25,  yO+128,  'Depth  (meters) ',  orient  =  90,  size  =  1.0,  color  =  0,  /device,  align  =  0.5 
skipgt:  print,  'skipping  plot' 

•  ■k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k-k'k'k'k-k'k'k'k-k'k'k'k'k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k'k'k-k'k-k'k-k'k-k'k'k'k 

r 

; COMPARE  RELATIVE  BATHYMETRY  MEASURE  TO  GROUND  TRUTH  BY  TAKING  RATIO  OF  RELATIVE  DEPTH  TO  TRUE  DEPTH 
; Ignore  areas  where  trudepth  =  0  (truedepth  data  don't  exist  at  those  points) 

truedepth (bad)  =  0 

good  =  where (truedepth  ne  0)  &  help,  good 

depth^ratio  =  fltarr ( samples ,  lines) 
depth  ratio  =  reldepth/truedepth 
depth^ratio (bad)  =  0 

rl  =  -0.5  ;min (depth_ratio (good) ) 
r2  =  1.0  ;max (depth_ratio (good) ) 

goto,  skip^rat  ;Not  really  useful 

window,  8,  xsize  =  samples,  ysize  =  lines,  title  =  ' (8)  Blue/Green  Band  Ratio  vs.  True  Bathymetry  -  '  + 
f  ilel9 
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im  =  bytscl(-l  *  depth  ratio,  min  =  rl,  max  =  r2) 
index  =  where ( land_gce_mask  eq  0) 
im (index)  =  255 
tv,  im,  order  =  1 

bar  =  bindgen(256) 

bar  =  reform(bar,  1,  256) 

bar  =  rebin (bar,  40,  256) 

xO  =  65  ;x0,  yO :  define  location  of  color  bar  within  the  image 

yO  =  30 

xl  =  xO  +  40;  +  1  ; xl ,  yl :  define  location  of  white  box  drawn  around  color  bar 

yl  =  yO  +  255;  +  1 
tv,  bar,  xO,  yO 

; Display  color  bar  on  the  image 

plo'  ,  bar,  /nodata,  /noerase,  color  =  0,  xrange  =  [0,  1],  yrange  =  [rl,  r2],  $ 
xticks  =  1,  yticks  =  2,  yminor  =  1,  charsize  =  1.0,  xstyle  =  1,  ystyle  =  1,  /device,  $ 
pos  =  [xO-l,  y0-l,  xl,  yl],  xtickname  =  ['  '],  yticklen  =  -0.05 

xyouts,  xl+25,  yO+128,  '  B/G  Band  Ratio  vs.  Depth  (m) ',  orient  =  90,  size  =  1.0,  color  =  0, 

align  =  0.5 

im  =  tvrd(true  =  1) 
im  =  reverse ( im, 1 ) 
im  =  reverse (im, 3) 
im  =  reverse ( im, 1 ) 

write_tiff,  outdir  +  band  ratio  vs  true  depth,  im 
skip_rat:  print,  'Skipping  ratio' 

•  -k'k'k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k 

r 

; REGRESS  RELDEPTH  VS  TRUEDEPTH  TO  FIND  ABSOLUTE  (DERIVED)  DEPTH  VALUE 

x  =  reldepth (good) 
y  =  truedepth (good) 

result  =  regress (x,  y,  SIGMA  =  sigmal,  CONST  =  const,  CHISQ  =  chisq,  yfit  =  yfit,  corr  =  corr) 
print,  'Equation  of  regression  line  (reldepth  vs.  truedepth) :  ' 

print,  'Regress  Result:  '  +  string (result)  +  'Regress  Constant:  '  +  string (const) 
print,  'Chi-squared:  ',  chisq,  '  Correlation:  ',  corr 


/ device. 
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YFIT  =  yfit2  ) 


print,  'Valid  regression?' 


degree  =  1 

Result2  =  POLY_F  (x,  y,  degree,  CHISQ  =  chisq2,  SIGMA  =  polysig, 
print,  'Equation  of  regression  line  (reldepth  vs.  truedepth) :  ' 

print,  'Polyfit  result:  '  +  string (result2 ) 
print,  'Chi-squared:  ',  chisq2 
print,  'Valid  regression?' 
help,  polysig,  yfit,  yfit2 

/Uncomment  to  see  the  scatter  plot  and  best-fit  regression  line  -  works,  just  slowly. 

window,  9,  xsize  =  900,  ysize  =  900,  title  =  ' (9)  Blue/Green  Band  Ratio  vs.  True  Depth  -  '  +  filel9 

x  =  reldepth (good) 

y  =  truedepth (good) 

good2  =  where (y  gt  -12) 
xx  =  x (good2 ) 
yy  =  y (good2 ) 


nele  =  n_elements (xx) 

index  =  indgen (nele/10,  /long) 

plot,  xx(index),  yy(index),  psym  =  3,  xtitle  =  'Blue/Green  Band  Ratio',  ytitle='True  Depth  (m) ',  title 
file,  $ 

yrange  =  [-15,  0],  xrange  =  [1.0,  1.2],  ymargin  =  [5,  4] 
oplot,  x (index),  yfit (index) 


slope  =  result2  (1) 
intercept  =  result2 (0) 


xyouts,  1.02,  -0.5, 
xyouts ,  1.02,  -1.0, 
xyouts,  1.02,  -1.5, 
xyouts ,  1.02,  -2.0, 


' Chi-squared 
' Correlation 
'Slope:  '  + 

' Intercept : 


:  '  +  string (chisq) 

:  '  +  string(corr,  format  =  "(f8.3)") 

string (slope,  format  =  "(f8.1)") 

'  +  string (intercept,  format  =  "(f8.1)") 


im 

im 

tv. 


=  tvrd ( ) 

=  255b  -  im 


im 


im 
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im  =  tvrd(true  =  1) 
im  =  reverse (im, 1) 
im  =  reverse (im, 3) 
im  =  reverse (im, 1) 

write_tiff,  outdir  +  regression_f ile,  im 

•  ■k'k-k'k-k'k'k'k'k'k-k'k-k'k'k'k'k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k-k'k-k'k-k'k-k'k'k'k-k'k-k'k-k'k-k'k-k'k'k'k-k'k 

r 

; USE  REGRESSION  RESULTS  TO  CORRECT  RELDEPTH  AND  CREATE  DERIVED  DEPTH 
absdepth  =  reldepth  *  result (0)  +  const 
rl  =  -18.0  ; 0 . 0 

r2  =  0.0  ; (-1 . 0) *min (absdepth) ; 0 . 0  ;2.0 

windov ,  10,  xsize  =  samples,  ysize  =  lines,  title  =  ' (10)  Derived  Bathymetry  -  '  +  filel9  /ABSOLUTE  IS 
THE  INCORRECT  TERM  —  DERIVED  IS  CORRECT 

im  =  bytscl (absdepth,  min  =  rl,  max  =  r2) 
index  =  where  ( land_gce__mask  eq  0) 
im (index)  =  255 

tv,  im,  order  =  1 

bar  =  bindgen(256) 

bar  =  reform (bar,  1,  256) 

bar  =  rebin (bar,  40,  256) 

x  0  =  65  +  50  ;x0,  yO :  define  location  of  color  bar  within  the  image 

yO  =  30 

xl  =  xO  +  40  ; +  1  ;xl,  yl :  define  location  of  white  box  drawn  around  color  bar 

yl  =  yO  +  255  ;6;  +  1 

tv,  bar,  xO,  yO 

/Display  color  bar  on  the  image 

plo-  ,  bar,  /nodata,  /noerase,  color  =  0,  xrange  =  [0,  1],  yrange  =  [rl,  r2],  $ 

xticks  =  1,  yticks  =  3,  yminor  =  1,  charsize  =  1.0,  xstyle  =  1,  ystyle  =  1,  /device,  $ 

pos  =  [xO-l,  y0-l,  xl,  yl]  ,  xtickname  =  ['  ,  yticklen  =  -0.05,  $ 

ytitle  =  'Blue/Green' 

xyouts,  xl+25,  yO+128,  'Derived  Depth  (meters) ',  orient  =  90,  size  =  1.0,  color  =  0,  /device,  align  = 

0.5 

im  =  tvrd(true  =  1) 
im  =  reverse (im, 1) 
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im  =  reverse (im, 3) 
im  =  reverse (im, 1) 

write_tiff,  outdir  +  f inal_depth_f ile,  im 

•  -k'k-k'k-k'k'k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k 

r 

x  =  absdepth (good) 
y  =  truedepth (good) 

result  =  regress (x,  y,  SIGMA  =  sigmal,  CONST  =  const,  CHISQ  =  chisq,  yfit  =  yfit,  corr  =  corr) 

print,  'Equation  of  regression  line  (absdepth  vs.  truedepth) :  ' 

print,  'Regress  Result:  '  +  string (result)  +  'Regress  Constant:  '  +  string (const) 
print,  'Chi-squared:  ',  chisq,  'Correlation:  ',  corr 

print,  'Valid  regression?' 

degree  =  1 

Result2  =  POLY_F  (  x,  y,  degree,  CHISQ  =  chisq2,  SIGMA  =  polysig  ,  YFIT  =  yfit2) 
print,  'Equation  of  regression  line  (absdepth  vs.  truedepth) :  ' 

print,  'Polyfit  result:  '  +  string (result2 ) 
print,  'Chi-squared:  ',  chisq2 

print,  'Valid  regression?' 
help,  polysig,  yfit,  yfit2 

/Uncomment  to  see  the  scatter  plot  and  best-fit  regression  line  -  works,  just  slowly. 

window,  11,  xsize  =  900,  ysize  =  900,  title  =  ' (11)  Derived  Bathymetry  vs.  True  Bathymetry' 

x  =  absdepth (good) 

y  =  truedepth (good) 

nele  =  n_elements (x) 

index  =  indgen (nele/10,  /long) 

plot,  x (index),  y (index),  psym  =  3,  xtitle  =  'Derived  Bathymetry  (m) ',  ytitle='True  Depth  (m) ',  title 
file,  $ 

yrange  =  [-15,  0],  xrange  =[-15,  0],  ymargin  =  [5,  4] 
oplot,  x (index),  yfit (index) 

slope  =  result2  (1) 
intercept  =  result2 (0) 

xyouts,  -14,  -1.0,  'Chi-squared:  '  +  string (chisq) 

xyouts,  -14,  -2,  'Correlation:  '  +  string (corr,  format  =  "(f8.3)") 
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im  =  tvrd ( ) 
im  =  255b  -  im 

tv,  im 


im  =  tvrd (true  =  1) 
im  =  reverse ( im, 1 ) 
im  =  reverse (im, 3) 
im  =  reverse (im, 1) 

write_tiff,  outdir  +  der_vs_true_scatter_plot_f ile, 

stop 

endfor 

end 


im 
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DERIVED  DEPTH  VS.  TRUE  DEPTH  PLOTS  AND  LINES  OF  BEST  FIT 


1010:  Cooslol  6c  Blue  Derived  vs.  True  Bothymetry 


1010:  Coo9tol  &  Green  Oerived  vs.  True  Bothymetry 


Derived  Bothymetry  (m) 


Derived  Bothymetry  (m) 


1010:  Blue  6c  Green  Derived  vs.  True  Bothymeti 


1010:  Green  6c  Yellow  Oerived  vs.  True  Bathymetry 


Derived  Bothymetry  (m) 


Derived  Bothymetry  (m) 


2010:  Coostol  &  Blue  Oerived  vs.  True  Bothymetry 


2010:  Coo9tol  Sc  Green  Derived  V3.  True  Bothymetry 


Derived  Bothymetry  (m) 


Derived  Bothymetry  (m) 


2010:  Blue  Sc  Green  Derived  vs.  True  Bothymetry 


2010:  Green  Sc  Yellow  Oerived  vs.  True  Bothymetry 


Cofftfotton:  0.2H 


5.5  -5.0  -4.5  -4.0 

Derived  Bothymetry  (m) 


Derived  Bothymetry  (m) 


3100:  Coostol  &  Blue  Oerived  v3.  True  Bothymetry 


3100:  Coo9tol  &  Green  Derived  vs.  True  Bothymetry 


Derived  Bothymetry  (m) 


Derived  Bothymetry  (m) 


3100:  Blue  Sc  Green  Derived  vs.  True  8othymetry 


3100:  Green  Sc  Yellow  Oerived  vs.  True  Bothymetry 


Derived  Bothymetry  (m) 


Derived  Bothymetry  (m) 


-■v  * 

■»*  - 

•  •  ’•  —v'-hl  -rv 

"  S 

4100:  Blue  &  Green  Derived  vs.  True  Bathymetry 
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TABLES  OF  EXAMPLE  CHI-SQUARED  AND  CORRELATION  VALUES 
FOR  EACH  BAND  COMBINATION 

1.  Coastal  &  Blue  (with  best  fits  highlighted) 


IMAGE 

ID 

meanSatEl 

Chi- 

Squared 

Correlation 

1010 

24.9 

4092798.6 

0.163 

1020 

26.5 

4081785.6 

0.169 

1030 

28.1 

4059720.3 

0.173 

1040 

29.8 

4043658.7 

0.19 

1050 

31.7 

3975405.5 

0.212 

1060 

33.7 

3878340.7 

0.254 

1080 

38.2 

3874253.7 

0.273 

1090 

40.8 

3740692.8 

0.319 

1100 

43.7 

3598962.8 

0.337 

2010 

46.7 

3520323.3 

0.372 

2020 

50.1 

3353757.3 

0.433 

2030 

53.6 

3377771.5 

0.425 

2040 

57.6 

3248184.9 

0.462 

2050 

61.6 

3216444 

0.472 

2060 

66.0 

3170698.1 

0.475 

2070 

70.2 

3187990.2 

0.48 

2080 

74.0 

3097967.8 

0.467 

2090 

76.9 

3107493.1 

0.483 

2100 

77.8 

3082345.9 

0.488 

3010 

76.3 

3105183.7 

0.489 

3020 

73.1 

3096221.6 

0.496 

3030 

69.0 

3119432.1 

0.493 

3040 

64.8 

3179182.1 

0.483 

3050 

60.6 

3215832.1 

0.475 

3060 

56.6 

3334628 

0.445 

3070 

52.9 

3375361 

0.434 

3080 

49.6 

3523816.2 

0.385 

3090 

46.5 

3585509.4 

0.359 

3100 

43.7 

3700396.6 

0.311 

4010 

41.1 

3768118.2 

0.28 

4020 

38.7 

3858040.5 

0.233 

4030 

36.5 

3893903 

0.212 

4040 

34.4 

3917829.2 

0.16 

4050 

32.5 

3875953.7 

0.139 
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4060 

30.7 

3757712.1 

0.127 

4070 

29.0 

3647051.5 

0.133 

4080 

27.4 

3520708.9 

0.113 

4090 

25.9 

3355284.3 

0.093 

4100 

24.5 

3020893 

0.038 

Coastal  &  Green  (with  best  fits  highlighted) 


IMAGE 

ID 

meanSatEl 

Chi- 

Squared 

Correlation 

1010 

24.9 

3911460 

0.262 

1020 

26.5 

3888263.3 

0.271 

1030 

28.1 

3882215.3 

0.277 

1040 

29.8 

3847912.4 

0.288 

1050 

31.7 

3755520.8 

0.313 

1060 

33.7 

3628560.2 

0.354 

1080 

38.2 

3546724.9 

0.39 

1090 

40.8 

3317653.4 

0.452 

1100 

43.7 

3132645.7 

0.475 

2010 

46.7 

3004689.7 

0.516 

2020 

50.1 

2762942.9 

0.573 

2030 

53.6 

2698408.7 

0.591 

2040 

57.6 

2585662.4 

0.611 

2050 

61.6 

2493514.7 

0.632 

2060 

66.0 

2406520.2 

0.644 

2070 

70.2 

2380175.1 

0.65 

2080 

74.0 

2331197.5 

0.651 

2090 

76.9 

2262424.3 

0.66 

2100 

77.8 

2242532.2 

0.665 

3010 

76.3 

2294378.5 

0.665 

3020 

73.1 

2291861.7 

0.663 

3030 

69.0 

2360830.5 

0.653 

3040 

64.8 

2429096 

0.636 

3050 

60.6 

2569034.7 

0.615 

3060 

56.6 

2718923.5 

0.588 

3070 

52.9 

2771285.8 

0.577 

3080 

49.6 

3030260.5 

0.519 

3090 

46.5 

3148438 

0.487 

3100 

43.7 

3326883.7 

0.434 

4010 

41.1 

3239760.9 

0.455 

90 


4020 

38.7 

3554008.1 

0.357 

4030 

36.5 

3580085.4 

0.345 

4040 

34.4 

3690842.1 

0.284 

4050 

32.5 

3740913 

0.23 

4060 

30.7 

3641605.2 

0.213 

4070 

29.0 

3543715.8 

0.212 

4080 

27.4 

3412602.3 

0.203 

4090 

25.9 

3204708.3 

0.143 

4100 

24.5 

2973309 

0.115 

Blue  &  Green  (with  best  fits  highlighted) 


IMAGE 

ID 

meanSatEl 

Chi- 

Squared 

Correlation 

1010 

24.9 

3837181.9 

0.294 

1020 

26.5 

3799491.4 

0.31 

1030 

28.1 

3826809.8 

0.295 

1040 

29.8 

3788817 

0.312 

1050 

31.7 

3674079.4 

0.34 

1060 

33.7 

3539442.7 

0.386 

1080 

38.2 

3449363.9 

0.419 

1090 

40.8 

3170962 

0.485 

1100 

43.7 

3013485.2 

0.512 

2010 

46.7 

2834575.7 

0.552 

2020 

50.1 

2629851 

0.601 

2030 

53.6 

2483369.3 

0.632 

2040 

57.6 

2372775.6 

0.652 

2050 

61.6 

2245770.1 

0.675 

2060 

66.0 

2157476.7 

0.689 

2070 

70.2 

2122878.3 

0.696 

2080 

74.0 

2099771.8 

0.698 

2090 

76.9 

2047377.9 

0.705 

2100 

77.8 

2019917.5 

0.707 

3010 

76.3 

2040766.2 

0.706 

3020 

73.1 

2064788.8 

0.703 

3030 

69.0 

2128616.8 

0.695 

3040 

64.8 

2209021 

0.682 

3050 

60.6 

2394974.9 

0.65 

3060 

56.6 

2537364.6 

0.626 

3070 

52.9 

2676606.9 

0.597 
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3080 

49.6 

2883049.7 

0.551 

3090 

46.5 

2994151.5 

0.521 

3100 

43.7 

3242079.3 

0.458 

4010 

41.1 

3074725.8 

0.499 

4020 

38.7 

3446563.2 

0.393 

4030 

36.5 

3455884.2 

0.391 

4040 

34.4 

3618164.1 

0.315 

4050 

32.5 

3676137.1 

0.268 

4060 

30.7 

3603561.2 

0.241 

4070 

29.0 

3515796.4 

0.23 

4080 

27.4 

3360508.4 

0.24 

4090 

25.9 

3203685.6 

0.232 

4100 

24.5 

2942699.1 

0.198 

Green  &  Yellow  (with  best  fits  highlighted) 


IMAGE 

ID 

meanSatEl 

Chi- 

Squared 

Correlation 

1010 

24.9 

4204121.8 

0.022 

1020 

26.5 

4199128.1 

0.056 

1030 

28.1 

4193524.9 

0.021 

1040 

29.8 

4195615.6 

0.008 

1050 

31.7 

4137792 

0.073 

1060 

33.7 

4088413.4 

0.121 

1080 

38.2 

4154217.8 

0.085 

1090 

40.8 

4059838.4 

0.161 

1100 

43.7 

4029843 

0.107 

2010 

46.7 

3892704 

0.212 

2020 

50.1 

3869278.8 

0.242 

2030 

53.6 

3803366.4 

0.278 

2040 

57.6 

3726302.1 

0.31 

2050 

61.6 

3732941.5 

0.309 

2060 

66.0 

3647327.9 

0.337 

2070 

70.2 

3558505.7 

0.374 

2080 

74.0 

3538108.6 

0.361 

2090 

76.9 

3526190.9 

0.36 

2100 

77.8 

3496716.8 

0.362 

3010 

76.3 

3537665.1 

0.371 

3020 

73.1 

3527937.5 

0.368 

3030 

69.0 

3792433.1 

0.267 
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3040 

64.8 

3710290.4 

0.322 

3050 

60.6 

4123412.3 

0.082 

3060 

56.6 

3981819.9 

0.21 

3070 

52.9 

3747042.4 

0.314 

3080 

49.6 

4140895 

0.012 

3090 

46.5 

4033670.6 

0.151 

3100 

43.7 

4077367.7 

0.07 

4010 

41.1 

3828814.7 

0.252 

4020 

38.7 

4043303 

0.095 

4030 

36.5 

4030102.2 

0.105 

4040 

34.4 

3891947.4 

0.18 

4050 

32.5 

3954476.1 

0.035 

4060 

30.7 

3821397.1 

0.045 

4070 

29.0 

3701158.5 

0.079 

4080 

27.4 

3565573 

0.068 

4090 

25.9 

3328975.3 

0.074 

4100 

24.5 

3061447.5 

0.064 
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